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Der kanonische Formalismus in entarteten Fällen 
Von Rudolf Haag in München 


Für Fälle, in denen die Gleichungen pa = O.L/d4a nicht nachden Geschwindigkeiten auflösbar sind, wird der 
kanonische Formalismus formuliert. Er,führt auf Systeme mit überzähligen Variablen. Ihre Reduktion und 
die Fragenach dem Sinn von Poisson klammerrelationen in solchen Systemen wird im 2. Teil besprochen. 


The “canonical formalism’’ is formulated for the case, that Ihe equations pa=9L]dg« cannot be solved for 
the velocities ga. It leads to a system with odd variables. T’he second part deals with the reduction of the system 
and discusses the meaning of the relations beiween the Poisson bracketts. 

Pour des cas, auxquels les &quations p, = OL/dq, me sont pas r6solubles vers les vElocites, le formalisme 


canonique est formul&. Il conduit ü& des systömes ü variables surnum£raires. Leur reduction et la question 
concernant le sens de relations en parenthese d’apres Poisson dans tels sysiömes est discutee dans la se- 
conde partie. 


Ha cıyyan, korna ypasHenua pa = 0L/d4a He PaspeImuMkL OTHOCHTENBHO CKOPOCTeÄ, dopMy- 
AuPyeTcH KAHOHNYecKnü hopMasınaMm. OH IPHBOAHT K CHCTEeMaM C HSÖBITOYHEIMH IIeEPeMeHHEIMM. 
Hx mpnsexenHue, a Takıke BONPOC O CMELIC.IE OTHOMIEHHH B CKODkaX IlyaccoHa B TAKHX CHcTeMax 
OÖCY:KAAMTCHA BO BTOPOH YacTH. 


1. Übergang zum kanonischen Formalismus 
Es sei die Bewegungsgleichung eines mechanischen Systems in Form des Variationsprinzips 
5 [L (9x, 9.) ds = 0; ER EN 


gegeben. Dann kann der Übergang zur Hamiltonfunktion und den kanonischen Be- 
wegungsgleichungen Schwierigkeiten bereiten, wenn die Legendre -Transformation 


oL 
= NEE (2a), 
H,= %I — L!) N a he per (2b) 


nicht gestattet, die Geschwindigkeiten g, als Funktionen der Impulse auszudrücken, wenn also die 
2 

Funktionaldeterminante 2 verschwindet. Gesetzt den Fall, sie haben den Rang (n—r), dann 
% 71B 

werden durch (2) gerade (n—r) der Größen 9, als Funktionen der p, und der restliche GeschwinA 

digkeitskomponenten ausdrückbar sein. Dafür werden r Beziehungen zwischen Impulsen und 

Ortskoordinaten bestehen: 

Z | BP 9) = 0; Er ER APahE RR EN) 


Wiebestimmt sich nun die Hamilton-Funktion in einem solchen Fall ? Zunächst bemerken 
wir, daß obgleich die 4, nicht vollständig durch 9,, q, ausdrückbar sind, der Übergang (2b) zu 
Ho(Px , 9.) stets möglich ist, denn das vollständige Differential von H,ergibt sich auch hierzu 


3 oL 
a a ir Pa ee te [AR 


_ Dawir.alle Größen als Funktionen von p,, 9, undetwag; (ü=1...r) ausdrücken können, H,sich 
"aber bei Konstanthalten der p,, g, nicht ändert, folgt, daß es unabhängig von 9; sein muß. D.h. 

aber: H,ist durch die p,, 9, allein ausdrückbar. Wegen der Beziehungen (3) ist dies allerdings 
“ auf mannigfache Weise möglich. Aus (4) folgen also nicht mehr die kanonischen Bewegungs- 


gleichungen. Vielmehr ist 


Er AH, E00, 
a -+ 4; = } . . . . . . . . Re RO ae I u 5 a ’ 
Ia OP OPx ) 
oH, ..0®, | 
ee a E65} 6) 5 
Pa ga A; 09x ( ) 


_ ) Wir verwenden durchweg die „Summenkonventionen“: über Indizes, die in einem Glied doppelt vor- 
mmen wird summijert. ; 
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wobei die A; zunächst unbestimmte Lagrangesche Multiplikatoren sind. Gehen wir von irgend 
einer der möglichen Darstellung für A, aus, so erscheinen die Gl. (5) in der üblichen Form ka- 


nonischer Gleichungen, wenn wir 
H= HD, ee (6) 


setzen und außer den zugehörigen kanonischen Gleichungen noch die Nebenbedingungen (3) 
fordern. Dieser formale Prozeß ist aber nur dann von Nutzen, wenn wir über die Multiplika- 


toren A; etwas aussagen können. Zu diesem Zweck beachten wir, daß nach (3) auch ®, ver- 
schwinden muß. Also 


Pe 


(HD, DD) 0 er (7). 


Soweit die Klammern {®,9,} nicht sämtlich verschwinden, ergeben sich hieraus Bestimmungs- 
gleichungen für die A,. Das System (7) wird aber im allgemeinen nicht sämtliche A; festlegen, 
sondern einige Gleichungen liefern, die von den Multiplikatoren frei sind. Sie sind entweder 
Identitäten oder neue Nebenbedingungen, etwa 


% (Dar Go) EN AdeT ra ee (8). 


Wiederholt man mit ihnen das Spiel, so kommt man schließlich zu dem vollständigen?) System 

von Nebenbedingungen und einer mehr oder weniger weitgehenden Festlegung der Multiplika- 

toren. 
Es gilt nun aber der folgende Satz: \ | 
Die aus(6)folgendenkanonischenGleichungenenthalten — zu- 

sammen mitdenNebenbedingungen(3) —alleAussagen über dieBe- 

wegung, welche in den Euler-Lagrangeschen Gleichungen des 

Variationsproblemsenthaltensind. 


Um dies einzusehen muß man sich nur davon überzeugen, daß aus den Gln. (5a) und (3) 
die 9, und A, als Funktionen der q,, 9, bestimmbar sind %). Dies zeigt man aber für die p,, indem 
man wieder durch Legendre-Transformation von H zu L zurückgeht. Man kommt dabei auf 
Gl. (2a), wenn man beachtet, daß die 9, alle unabhängig sind, sofern nur die Gln. (5a) und (3) 
berücksichtigt werden und nicht auch die Bewegungsgleichungen (5b) und die daraus fließenden 
Kenntnisse über die A; (Gl. (7) usw.). Daß schließlich die A; durch die 9, ausdrückbar sind, er- 
kennt man am einfachsten, wenn man die Nebenbedingungen (3) in einer, nach r Impulsen 9,, 
aufgelösten Form einsetzt: 


B—kP,ga)=0 i=l...nr, s=er+l...n 
H, = Ho(Pe, 4%): 


TEE ER ne a rep 


2 Dann liefert (5a) direkt: 
RI; 


Wenn nun die Bedingungen ®; = 0, ®,=0 usw. noch einige der Multiplikatoren unbestimmt 
lassen, so folgt, daß jede Wahl dieser Multiplikatoren mit den Bewegungsgleichungen verträglich 
ist. DieBewegungsgleichungenplusAnfangsbedingungenbestim- 
mendieFunktionenq(s)nicht vollständig. Einige bleiben willkürlich wähl- 
bar. Wenn andererseits der physikalische Sachverhalt ein wohlbestimmter ist, so bedeutet dies, 
daß die Größen g und snichtsämtlicheine selbständige physika- 
lische Bedeutunghaben. 

Beer, Ein häufiges Beispiel ist der Fall der „Parameterinvarianz‘“: Das Extremal- 
Br prinzip ist invariant gegen eine beliebige Änderung des Zeitparameters s—f(s). In diesem 
Br Fall ist s physikalisch belanglos. Dies wirkt sich auf die Hamiltonfunktion so aus, daß sie 
nur bis auf einen beliebigen Faktor bestimmt ist. Ein anderes Beispiel: | 


ı r 
L= = El RP ET EDER Pa en | (9). 
‘ Hier besteht Invarianz gegen die Transformation 


G>u- OR s>3=fß). . 


2) Unter {4, B} verstehen wir die Poissonklammer von A und B. 
.08A0B 0Ao0B 


®) Das System ist vollständig, wenn alle weiteren Differentiationen zu Identitäten f 
*) Denn dannliefert (5b) nach Elimination von px und A; die Euler-La gra i g ruken Glasur 


ee 
, nd pr 
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In diesem Fall haben also die Verhältnisse der Koordinaten und das Differential ER invariante 
Bedeutung. (4)? 


Esist vielleicht von Nutzen, die obige, etwas formale Schlußweise noch zu ergänzen. Wir 
beschränken uns auf den Fallr=1. Es soll also die Funktionaldeterminante 


®L 

09% Ge 

den Rang (n— 1) besitzen. Anders ausgedrückt: Es existieren (bis auf einen gemeinsamen 
Faktor bestimmte) Funktionen c, (9, 9) so daß 


Opa Opa 
Ser = ee ie DE Et LE, er ee 10 
. 7 09% 098 
Der Differentialoperator 
(] 
EX = Cx 24, mer, Do Mae Er Marie ee a el are (11) 


hat dann die Eigenschaft, daß er alle Größen vernichtet, diesich als Funktionen von Pa, 9. allein 
ausdrücken lassen. Die c, sind einerseits die Unterdeterminanten einer Zeile der Funktional- 
matrix, andererseits stehen sie in engem Zusammenhang mit der Funktion ® von Gl. (3) 
= ar I Bewegungsgleichungen 
DeLe 04: oL 
20,0, 4% ee 
erhält man durch Multiplikation mit c, eine Differentialgleichung 1. Ordnung: 
| öl oL 
0x 090g Ie — © gr — 
Diese Gleichung ist identisch mit ®=0. Da nämlich 
D % ’ 1) —=0 
identisch in g, und q, erfüllt ist, folgt durch Differentiation nach gg: 


und durch Differentiation nach 95 En 
; 0D 09x 

N EL er te (15). 

Op« 2gp " 2gp 

Vergleich von (14) und (10) zeigt zunächst, daß bis auf einen willkürlichen Proportionalitäts- 

faktor | 


= N ER EN (16) 
Opa ” 

Setzt man dies in (13) ein und beachtet (15) und die Bewegungsgleichung En — p,„ so entsteht: 
& 


od. od . . 
— + — Hp =d=N. 
EEE + a 
Aus dieser Bedeutung der Gl. (13) ersieht man insbesondere, daß die linke Seite eine Funktion 
von p, und q, allein ist, nämlich die Größe — {H,, ©}. 
Durch Anwendung des Operators X läßt sich dies auch direkt verifizieren. Esist 
ö @L_, RER AR @oL @L\ 
DE CHE kr 7 ki (areas 10) + er0e (age dar — rn) 


=R > dB Zen) : 
2 04, \0gu 098 99x 
Der zweite Term verschwindet wegen der Antisymmetrie der Klammer. Den ersten Term schreiben wir 
a N #017 ..06% 
2 nr (or Bu i) a P dge 

Beide Ausdrücke verschwinden wegen (10). Schließlich verschwindet der dritte Term, wenn man 
—= Ca(P, q) setzt, wegen der Eigenschaft des Operators X. 

| ” | 14* 


4 


durch (19a) eingeschränkt. 


“ De Su Pr Dr 
a it 7° ee Re LS, 

EEE 8ER j< * 
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Wir haben also den folgenden Sachverhalt, wenn wir q,, 9, als unabhängige Variable be- 
trachten: Die Gleichung ® = 0 ist identisch erfüllt. ® = 0 ist gleichbedeutend mit der Dif- 
ferentialgleichung 1. Ordnung (13). Die Bewegungsgleichungen (12) lassen Auflösung nach 


= I (9a, Ir H); (=2...n) 
zu An Stelle einer Aussage über g, hat man zunächst die Gl. (13). Ihre Differentiation gibt 
wiederum keine Aussage über {,, sondern eine weitere Gleichung 1. Ordnung, da die linke Seite 
als Funktion von p,, q, darstellbar ist, und 9,= = nur die ersten Ableitungen enthält. 
& 

Esgibtnun zwei Möglichkeiten : Entweder ® ist ebenfallsals Funktionder pund g darstellbar. 
Dann liefern die Differentiation wieder eine Gleichung 1. Ordnung. Oder aber ® enthält neben 
p,q noch q,. Dann liefert die Differentiation die gesuchte Aussage über g,. Bei der fortgesetzten 
Differentiation der Gl. ® = 0 ergibt sich also eine Folge von Differentialgleichungen 1. Ordnung, 
d.h. Einschränkungen für die Anfangsbedingungen. Die Folge endet entweder dadurch, daß 
schließlich eine Identität resultiert. In diesem Fall ist die Funktion 9, willkürlich wählbar. Oder 
es ergibt sich schließlich eine Gleichung 2. Ordnung durch die q, als Funktion der g,, 9, fest- 
gelegt wird. Im kanonischen Formalismus ist dies der Fall, in dem der Multiplikator bestimm- 
bar ist. 

Wir wollen diese Beziehungen nochmäls an Hand des Beispiels (9) verfolgen. Aus 


L= 5 [192 — (44) 
folgt 

Pa = YuIa — (UuIu) 9a - 
Zwischen den 7,,g. besteht die Beziehung 


Legendre-Transformation führt zu der Hamiltonfunktion 


2 1» 
| B=mu—ı=L-gn 
Die vollständige Hamiltonfunktion lautet also 
_1Ps 
a 2q + y a (Pa 1) Ta Ru Zei Eu (18). 


Differentiation von (17) liefert 


oder 


a von (19) liefert nichts Neues mehr. Der Multiplikator bleibt also un- 


Gehen wir andererseits von den Bewegungsgleichungen aus: 
Qudx u (gu du) 1x 7 2(qu9u) da HW7 2(9) I = Dr ehe EM (20), 


so folgt durch Multiplikation mit 3 = 9. die Gleichung 1. Ordnung 
& 


N) DE ALLEN (19a). 
Sie entspricht, wie wir oben allgemein zeigten der Gl.® = 0 (19). Differentiation der Gl. (19a) 


liefert nichts Neues, i N 5 
al ehreiben eues, sondern eine Konsequenz von Gl. (20). Die Bewegungsgleichungen lassen 


re 2(Q, 0 
= (q) + A 1 RR a a LEN EA + (20a) 


wobei A eine willkürlich wählbare Funktion ist. Die möglichen Anfangsbedingungen werden 


r e. 
5. 
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2. Reduktion der Variablensätze und Poisson klammerrelationen 
Gegeben eine Hamiltonfunktion H(p,, 9.) und eine Anzahl von Nebenbedingungen 


€ | Pier 1) 0. 
Wir wollen übergehen zu einer geringeren Zahl von Variablen die unabhängig sind °). 
Es bestehe also eine umkehrbare eindeutige Beziehung 


Pa: fa *> P,, 0, ®;. 
H(px; > Zur. Q,)- 


Die Bedingung dafür, daß die kanonischen Gleichungen in den beiden Variablensystemen 
äquivalent sind, ist, wie man leicht nachrechnet: 


Ihr zufolge entspricht 


oP, OP (r) od, 
FF + 0, EN ar rn (1a), 
op, 04x __.4@) oB, 
ae ee (1b). 


ae r se mit Multiplikatoren 1, welche ebenso 
& & 
wie die A} willkürlich sind. Die Gleichungen sind ungeschickt, da in dem ersten Glied links je- 
weils die alten, im zweiten Glied die neuen Größen als unabhängige Variable fungieren. 

Man kann sie umformen ®) in 


BREI Bl 0; Wa ee a Ey (2), 


wobei die „Lagrangeschen Klammerausdrücke“ definiert sind als 


Dazu kommen entsprechende Gleichungen für 


OP ga 09x OP 
[4, Bl = I BER N (3). 


In diesenGleichungen wird nur nach den neuen Koordinaten differenziert. Meist ist die Benutzung 
der alten Größen als unabhängige Variable bequemer. Hier ergibt sich das Gleichungssystem ”), 2) 


I) 9, = bs 7 Dun C;; 

{P,, P}= A Oi; 

N (4) 
{P,8)}=190,, 

0, 0}= w0C;; 


{®; P D} = ER RE a EFT REIT EE e, (5). 
9, 0, 


Kombiniert man die Gl. (1) mit den entsprechenden Gln. für REPE so erhält man die Be- 
& %&% 
dingungen ausgedrückt in den Poissonklammern hinsichtlich P, @ und in den Lagrange- 
Klammern hinsichtlich p, q. Sie sind jedoch weniger bequem, so daß wir sie hier nicht an- 
schreiben wollen. Im Folgenden knüpfen wir an (4) an. 
Wenn überhaupt eine Reduktion der Variablen möglich ist, also eineLösung für P, Q, A, 
win (4) existiert, so gibt es auch eine mit A, u = 0. Wir brauchen dazu von der als gegeben an- 


genommenen Lösung von (4) nur überzugehen zu 

: | PR,=B—-% 9; 9,=9,— u ©;. 
Dies ist von Bedeutung für die Beurteilung von Poissonklammerrelationen im System der 9, 
9a. Bezeichnen wir mit {4, B}, die Poissonklammer der Größen A und B im System der Pa» Ia 
und mit {A, B},; dieKlammer im System P,, Q,, so hat{4, B}; einen eindeutigen physikalischen 
Sinn, dagegen {AB}, nicht, denn dieselben physikalischen Größen A, Bkann man zufolge der 


5) Ist die Anzahl der ©; ungerade, so kann dies natürlich nicht vollständig gelingen. Indessen gelten 
die folgenden Überlegungen auch bei einer teilweisen Elimination der Nebenbedingungen. 


mit der Abkürzung 


6) Indem man (1a) mit a (1b) mit — an multipliziert und addiert ergibtsich z. B.die1. Gl. (2). 
?) Z.B. Multiplikation von (la) mit — Ze, von (1b) mit Fr und Addition. 
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Nebenbedingungen ©; = 0 in mannigfacher Weise als Funktion der 9, 9x darstellen. Physi- 
kalisch sinnvoll ist eine Poissonklammer {AB}, offenbar dann, wenn A und B in solcher 
Weise als Funktionen der p,, 9. dargestellt sind, daß 


{4, Bi, = {4, Bir: 


Dies ist aber nach dem obigen gerade dann der Fall, wenn A und Bals Funktionen der P, Q 
allein dargestelltsind. Diese sind jedoch dadurch ausgezeichnet, daß ihre Poissonklammern mit 
den ®, verschwinden. Dasselbe müssen wir daher für A und B fordern und wir haben den Satz: 
Damitdie Poissonklammerrelation {A, B};=(C physikalisch sinnvoll ist, müssen AundB 
in einer solchen Form angeschrieben werden, daß ihre Poissonklammern mit den Neben- 
bedingungen verschwinden: 
{A, B;h = {B, ®;h hr 
Eingegangen am 8. Sept. 1951. 


Fehlerabschätzung bei der ersten Randwertaufgabe 
bei elliptischen Differentialgleichungen 


Von ZL. Collatz in Hannover 


Der Satz vom Randmaximum für subharmonische Funktionen, der unter gewissen Voraussetzungen auch 
für allgemeinere lineare elliptische Differentialgleichungen gilt, dient bei der 1. Randwertaufgabe zur Fehler- 
abschätzung für eine Näherungsfunktion, sowohl wenn diese wie z.B. beim Trefftzschen Verfahren nur die 
Differentialgleichung erfüllt, als auch, wenn sie wie z.B. beim Ritzschen Verfahren nur der Randbedingung 
genügt. Die Abschätzungen werden angewendet auf ein Verfahren, das auf eine Idee von R.v. Mises zurück- 
geht, bei welchem man Differentialgleichung und Randbedingungen simultan annähert. 

For the first boundary value problem, the boundary maximum theorem for subharmonie functions, that is, 
under certain conditions, valid also for more general linear differential equations of the elliptic type, is used 
for the estimation of the error of the approwimate function. T’he approximation may only satisfy either the 
differential equation (Trefftz’method) or the boundary conditions ( Ritz’method). The estimations are 
applied to a method that — following an idea of R.v. Mises — approzimates the differential equation and 
the boundary conditions simulianeously. n, 

Le theoröme du maximum marginal pour des fonctions sous-harmoniques; lequel, ü certaines supposi- 
tions, s’applique aussi a des Equations differentielles lineaires, elliptiques plus gen£rales, au premier pro- 
bleme de valeur marginale, sert & l’estimation des erreurs pour une fonction approximative, et au cas que 
celle-ci, comme p.e.a la methode de Treffz, ne correspond qu’& l’&quation differentielle, et au cas que, comme 
p.e. a la methode de Ritz, elle ne satisfait qu’& la condition marginale. Les estimations sont appliquses & 
une methode, qui a son origine dans une idee de R. v. Mises, ülaquelle on approche l’&quation differen- 
tielle et les conditions marginales simultan&ment. 

Teopema EpaeBoro MakCHMyMa 714 CyÖTapMOHHYecKHX PyHKUHK, MeÄcTBHTeNsHaA IIPH 
U3BeCTHBIX IIPEMIIOCHLIKAX TAKIKE IA 6oNee OÖLIHX IHHEHHHX SILIHIIHYeCKHX AubhepeHuaAsHEIX 
yPaBHeHHü, CIIyKUT B IIePBOH KpaeBofi 3anaye IA OMeHKH OMMHÖOR MIA IIPHÖNFKEeHHOH DYHK- 
IHN; 9TO UMeeT MecTO KAK B TOM CIY4Ae, KOT TOCHENHAA — KAK HAIP. IIPH CIHOCoÖe u 
— B5INONHAeT TOANBKO AydhepeHmmanbHoe yPpaBHeHHe, MaK U B TOM, KOTI& OHA — Kak HAIp. 
ıpu cmocode Purua — yN1oBIIeTBOpAeT TOIBKO KPaeBOMy YCAIOBHW. ÜÖHeHEH IIPHMeHADTCH K 
cITOCoÖy. OCHOBBIBAWINeMycA Ha unee P.d&. Museca pm KOTOPoM AudhepermmansHoe YbaBHeHHe 


Hu KpaeBhle YCIOBHA COBMECTHO IIPHÖNMKAMTCA. 

In Nr. 1 bis Nr. 3 wird zunächst eine allgemeine Fehlerabschätzung für die 1. Randwert- 
aufgabe bei elliptischen Differentialgleichungen auf Grund des Satzes vom Randmaximum auf- 
gestellt, in Nr. 4. das Trefftzsche Verfahren!) etwas ausführlicher behandelt; in Nr. 6 bis Nr.9 
wird eine auf R. v. Mises zurückgehende Methode der gleichzeitigen Annäherung von Differen- 
tialgleichung und Randbedingungen dargestellt und auf sie die Fehlerabschätzungen von Nr. 2 
und Nr. 3 angewendet. Diese Methode ist auch bei allgemeineren Randwertaufgaben anwendbar, 
und daher sind Nr.6 bis Nr. 9 so ausgeführt, daß sie (abgesehen natürlich von den Verweisen auf 
den ersten Teil) auch für sich ohne Lesen von Nr. 1 bis 5 verständlich sind. Deshalb ist auch in 
Nr. 6 bis Nr. 9. von der in Nr. 1 bis 5 verwendeten abgekürzten Schreibweise des Fortlassens der 
Summenzeichen (die sich bei Beweisen gut bewährt) kein Gebrauch gemacht. 


1. Bezeichnungen 
Vorgelegt sei die lineare Differentialgleichung zweiter Ordnung 
Lluler Sn RE N) 


; n n 
Dit L[u]= = Paul Up 3 bu;tew. 


1) E. Trefftz, Ein G - h x 
NE Eh er H ges zum Ritzschen Verfahren, Verh. 2. internationaler Kongreß f. tech. 


Rn lt efftz, Konvergenz und Fehlerabschätzung beim Ritz schen Verfahren, Math. Annalen 100 (1928), 


£ 
ü 
{ 
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Dabei ist u(2, ..., 2,) eine gesuchte Funktion der unabhängigen reellen Variablen &,,..., ©, 
in einem Gebiet B-+-I/' des n-dimensionalen Raumes, die Gjp d;, c, rsind in B-++-I' gegebene, etwa 
stetige Funktionen derselben Argumente, und in Gl. (2) wie auch im folgenden sind (wie in der 
Tensorrechnung üblich) die Summenzeichen fortgelassen, indem in einem Produkt über einen 
zweimal auftretenden lateinischen Index von 1 bis n zu summieren ist. Bei@und b dienen an- 
gehängte Indizes 7, k nur zur Unterscheidung von Funktionen, bei %, v, w, 2 bedeuten sie partielle 
Ableitungen: 
BR ou Beru 
> : 0%; £ Die” 0%; 0% ; 


Die Differentialgleichung sei elliptisch, d.h. die Matrix der G;, (die als symmetrisch vorausgesetzt 
werden kann: q,.— a,,) sei positiv definitin B-H7'. Fernerseic>0. ®seiein offenes beschränk- 
tes Gebiet und ®-+-/’ das durch Hinzunahme der Randpunkte aus ® hervorgehende abgeschlos- 
sene Gebiet. /’sei dabei eine abgeschlossene zusammenhängende (n—1)-dimensionale stückweise 
- glatte?) Hyperfläche. e 
Es gilt der Gaußsche Integralsatz für eine in 8-7’ stetige, mit stetiger Ableitung nach 

einer Variablen ; versehene Funktion «(2,,..., 2,): = 


[3a = [eos 6 mar En ae Be Pe ee 
8 ;; 


0%; 


Dabei bedeuten dr=dx,dx,.... dx, das Volumenelement, df das Oberflächenelement auf T' und 
v die innere Normale. Die Normale darf dabei auf gewissen (n—2)-dimensionalen Mannigfaltig- 
keiten nicht definiert sein. Fe 34 

Nun seien auf I’ die Randwerte % von u vorgeschrieben, wir nehmen der Einfachheit halber 


u auf I’ als stetig an (auch bei anderen Funktionen v,.... soll ®, .... jeweils die Randwerte auf 
T’ bezeichnen). Die 1. Randwertaufgabe verlangt dann die Bestimmung von v aus (1) und 
5 In ET RS RE TERN NE 


‘ 


| 2. Ein Hilfssatz 
_ InNr.3 wird ein Hilfssatz benötigt, der in ähnlicher Form bei Courant-Hilbert,a.a.O., 
S.275 und bei Bateman’°)steht; da er aber hier etwas allgemeiner gebraucht wird, sei sein Beweis 
kurz angegeben. i 
Hilfssatz: Für eine nichtkonstante, mit stetigen partiellen Ableitungen bis zur zweiten Ord- 


_ die Hyperkugel X, mit P, als Mittelpunkt und R=P,P, als Radius gilt dann: In 


dem Öberflächenpunkt P, ist w(P,)=M, sonst ist im Innern und auf der Oberfläche 
der Hyperkugel w<M. Um den Punkt P, wird eine Hyperkugel X, von einem Radius k<R 


- gelegt. Nun werde vorübergehend der Punkt P, als Anfangspunkt des SR %,-Koordinaten- 
systems gewählt und die Entfernung von ihm mit r bezeichnet, r=%, 2. Die Hilfsfunktion 
2 g=wHö [etr—e-*#’], wo aund ö positive Konstanten seien, ist nun — M im Mittelpunkt P, von 
K,und <M auf der ganzen Oberfläche ® von K, bei passend gewähltem o, 6; denn auf dem Teil | 
yon d, der außerhalb X; liegt, ist die eckige Klammer negativ und auf dem mit X, gemeinsamen 
- — Teil von Bist w<M, und bei gegebenem « kann man 650 klein wählen, daß dann dort auch 2<M 
bleibt. Also: znimmt im Mittelpunkt von K, den-Wert. M ‘an und hat auf der ganzen Oberfläche 


N 2 | Fri Courant -D. Hilbert, Methoden der mathematischen Physik, Bd. 2, Bln. 1937, 8. 228. 


a ®) H. Batema n, Partial differential equations of mathematical physics, New York 1944, S. 135. 


u » + 
r% 3 “r je Pi 
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von K, Werte z<M, es muß also z mindestens in einem Punkt P; innerhalb von K, ein differen- 
zierbares Maximum annehmen. 
Nun rechnet man leicht aus: 


Liz) = L[w]+6 L[e"']—c de *" 

mit er {L[e"’] —ce"R'} = — 40? 0,,%;%4420(0;;+b; 2) +c— cert’—-M, 

Hier hat «? einen negativen Faktor, und das letzte Glied ist wegen c>0 nichtpositiv; durch 
passende Wahl eines genügend großen & kann daher erreicht werden, daß die geschweifte Klammer 
an der Stelle P, negativ (Null ausgeschlossen) wird. Wegen Z[w] <0 muß also L[z] <O sein 
an der Stelle P,. Nun kann aber L[z] an der Stelle P, unmittelbar angegeben werden. Dort 
hat z ein differenzierbares Maximum, es ist ,—0 (für j=1,2,:. „.N), und die negativ ge- 
nommenen zweiten Ableitungen bilden eine positiv definite oder semidefinite Matrix 


| — 2466,20 
für beliebige reelle &,. Es wird also 
L[2]p, =—% 26 +6 2. 
Da sowohl die a,; als auch die —2;; die Koeffizienten positiv (semi-)definiter quadratischer For- 
men sind, ist nach einer Formel von Fejer®) 
— Ag 2jk > 0. 

Da z(P,)=M ist, muß das Maximum von z, also z(P,) ZM sein; nach den Voraussetzungen ist 
cz=0 im Falle c=0 und cz(P,)=0 im Falle c>0, da dann M 20 vorausgesetzt war; mit- 
hin ist Z[z]r, >0, während die vorangehende Betrachtung Z[z]r, <0 ergab, womit ein Wider- 
spruch vorliegt. 3 

Corollar: Durch Übergang von w zu —w sieht man, daß der Hilfssatz bestehen bleibt, wenn 
man bei L[w] und M die Ungleichheitszeichen zZ, < vertauscht und als M das Minimum von w 
nimmt. 


3. Allgemeine Fehlerabschätzung 


Aus dem Hilfssatz ergibt sich sofort eine Fehlerabschätzung für Näherungslösungen v der 
ersten Randwertaufgabe (1), (4). Dabei setzen wir voraus, daß v von der Randwertaufgabe (1), 
(4) wenigstens eine der beiden Gleichungen (1), (4) erfüllt und unterscheiden somit zwei Fälle. 

1. Fall: v erfüllt die Differentialgleichung ZL[w]=r exakt, braucht aber nicht die vor- 
gegebenen Randwerte % anzunehmen. Die Fehlerfunktion w=v—u genügt dann der Rand- 
wertaufgabe L[w]=0, w=v—u, wobei die Randwerte ö mit gegebener Näherungslösung 
ebenfalls bekannt sind. Mit Z[w]=0 sind beide im Hilfssatz und im Corollar auftretenden 
Voraussetzungen Z[w]<0 und Z[w] >20 erfüllt. Wenn w überhaupt nichtnegative Werte 
annimmt, ist das Maximum M von w also M=0, und nach dem Hilfssatz nimmt w den 
Wert M nur auf dem Rande I’ an; es ist also dann War Wmaz 20. Nimmt w überhaupt 
nichtpositive Werte an, so ist entsprechend nach dem Corollar win = Umin<0. Nimmt 


‚also w den Wert O an, so ist Oin SW < Üpaz, und es ist Yin SA mas 0; nimmt w nicht den 


Wert 0 an, hat also w als stetige Funktion ein festes Vorzeichen, z.B. w>0, so ist O< w< Wax: 
Wir können die verschiedenen Fälle auch nach dem Verhalten von ®ö an Stelle von w ord- 
nen, was für die Anwendungen bequemer ist, und dann ergibt sich die x 
Fehlerabschätzung: Bei der Randwertaufgabe (1), (2), (4) mit c 20, a,, positiv definit 
sei v eine die Differentialgleichung L[v]=r erfüllende Näherungslösung. Von der Fehlerfunktion 
w= v—u sind die Randwerte w=®—U bestimmbar. Nimmt ö den Wert Null an, so gilt in B die 


Fehlerabschätzung > RER ) 
min Mate Er A Ne EP Sr g 


Nimmt w nicht den Wert Null an, und hat also wein festes Vorzeichen, so hat auch die Fehlerfunk- 


tion w in B + I’ überall dasselbe Vorzeichen, und für den Betrag von w gilt 


ol & Ta + E E (6). 


Im Falle c=0 gilt (5) unabhängig davon, ob w den Wert Null annimmt oder nicht. 

Vom Standpunkt der praktischen Analysis aus liegt hier ein günstiger Fall für eine Fehler- 
abschätzung vor. Solche Abschätzungen fallen bei anderen Differentialgleichungsaufgaben oft 
umständlicher und ungenauer aus, während die Abschätzungen (5) und (6), sofern man sie als 
für den ganzen Bereich ® gültig auffaßt, die bestmöglichen sind; der Fehler erreicht ja tatsäch- 
lich die angegebenen Schranken auf dem Rande I" 

2. Fall: Die Näherungslösung v erfüllt die Randbedingungen (4) exakt, braucht aber nicht 
der Differentialgleichung (1) zu genügen. Dieser Fall tritt z.B. auf, wenn man das Ritzsche Ver- 


4) L. Fej6r, Über die Eindeutigkeit der Lösung der li j i ialglei ‚zwei 
Ordnung, Math. Z.1 (1918) 8. ee g der linearen partiellen euren zweiter 
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fahren durchgeführt hat; er besitzt daher besonderes Interesse. Es wird jetzt eine Fehler- 
funktion Lfvy]— L[u]=L[v]—r auftreten; nun verwende man zwei Hilfsfunktionen 9» 9a 


derart, daß 
Em Katze rw Lioitliojer..... 2.0.0.5 (7) 


in ganz ®+-!' gilt. Man stelle etwa Z[v] auf, bilde mit beliebigen Funktionen qi und g2 die 
Ausdrücke L[gi] und Z[g5] und suche durch Kombination von v mit ,—=c,gf bzw. P=nR 
die Ungleichungen in 8+-/’ zu erfüllen. Stets lassen sich derartige Funktionen g,, 9, angeben, 
wenn man eine Funktion z mit Lg] <—A<0 in B kennt; z.B. für den Potentialausdruck 


N 
ou AN, 
L[u])=— = dal kann man 2= 5, 2 x; verwenden; man braucht dann nur in ;=t-+Y;2 
= I= 


die Konstanten y, und y, passend zu bestimmen. Dann ist der Hilfssatz auf die Funktionen 
w=%+91— u, bzw. v+ 9,— u anwendbar. Die Werte der Funktion v+g,— u sind am Rande be- 
kannt, der Maximalwert am Rande sei M,, entsprechend sei M, der Minimalwert von v+,—u 
am Rande. Ferner setzen wir 
el, im Falle c=0 5 
#1] Max(M,, 0) im Falle c>0 
nr j&M, im Falle c=0 
re 0) im Falle c>0. 
Dann besagt der Hilfssatz, bzw. das Corollar 


vH N—UuSsm Vo —USh 
und damit erhalten wir die Schranken für «: 
OMU SUSI TUE ee re (8). 


4. Aufstellung einer Näherungslösung v und Trefftz sches Verfahren 


Von den zahlreichen Möglichkeiten, eine Näherungslösung ® aufzustellen, seien hier nur - 
die folgenden genannt?). 
Für v wird der Ansatz gemacht 


BI U U Fe a Ri ee og (9). 

Dabei möge über einen in einem Produkt zweimal auftretenden griechischen Index 0, 0,... von 

1 bis p summiert werden. Es sind v,, v1, .. . ‚u, fest gewählte Ansatzfunktionen, welche die in- 
homogene, bzw. die homogene Differentialgleichung (1) erfüllen: 

Livl=r, L[%]=0 Bern). 2r(10): 


Dann erfüllt in (9) vfür beliebige Werte der Parameter a, die Differentialgleichung (1). 
Ein naheliegender Weg zur Bestimmung der a, ist die Fehlerquadratmethode; dabei soll 
v am Rande möglichst gut die gegebenen Randwerte u annähern, d.h. es soll 


El Mina a en (11) 


einen Kleinstwert annehmen. Dann sind 
10J 2 ER. 
— — = | (dt % d% — u) d, df=0 (o=1, 39%;9) 712) 
2.00, 
r 
plineare Gleichungen für die Unbekannten a,. Für die zugehörige Funktion v.nach (9) kann man 
dann nach Nr. 3 eine Fehlerabschätzung durchführen. 
Trefftz beschritt für die Potentialgleichung einen anderen Weg, bei welchem man zu- 
gleich für ein bei den Anwendungen öfter gebrauchtes Integral Schranken erhält. 
Die Differentialgleichung wird jetzt als selbstadjungiert vorausgesetzt, d.h. es sei 
00,7 
b b; 2% 0% 
Wie beim Ritzschen Verfahren spielt auch beim Trefftzschen Verfahren ein Integral J 
(eine Verallgemeinerung des Dirichletschen Integrals) eine wichtige Rolle: 


- Ilv, o]= [uw tow)a=J[wv]) 2.0 .: (13). 

3 « 
Es wird J[v, v]=J[v] gesetzt. Wendet man den Gaußschen Integralsatz (3) auf die Funktion 
a= VA, Wr, an, SO erhält man eine der Greenschen Formeln, wobei man in der entstehenden 
3 5) Eine Aufzählung verschiedener Arten findet sich z.B. Z. angew. Math.Mech.,Bd.31(1951), 8.230— 236. 


(Summation über %). 


Ai 


& 
E 
5 

5 


x 
2 
hr 
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Formel die rechte Seite bekanntlich durch Einführung der Konormalen®) o, die mit einer Größe A 


zusammen durch 
4;, c0s(v, ©) A cos (0, %,) (Ein) 


bestimmt wird, umformen kann: 


0) 
[ vrtwIae= 415 w]-+ fra ee. (14). 
8 R 
Zur Bestimmung der Parameter a, im Ansatz (9) wird jetzt die Forderung aufgestellt, daß der 
Fehler w=v—u in der durch das Integral J vermittelten Metrik gemessen zum Minimum ge- 
macht werden soll, also 
J[w] =Jfo—u]=Min. .. vr» Fan (15). 
Die Forderung (15) ergibt bei: dem Ansatz (9) 

10J[o+% %— u] 

2 da, #® 
Nach (14) folgen 


=J[o+ %W—u, u, ]=J[w, v,]=0 Ver li 


«Au z a Were 


r F (über og summieren) 


als p lineare Bestimmungsgleichungen für die a, (die Größen v,, %,, u sind auf /’ bekannt). Hat 
man aus diesem Gleichungssystem die a, berechnet, so kann man für die mit diesen a, nach (9) 
gebildete Funktion v wieder nach Nr. 3 eine Fehlerabschätzung durchführen. 


5. Beispiele 
Die beiden folgenden Beispiele dienen lediglich zur Erläuterung der Methoden, es ist dabei 
nicht. auf große Genauigkeit Wert gelegt. 
I. Als einfaches Beispiel werde die stationäre Temperaturverteilung «(z, y, z) im Innern 
eines Würfels |x]<1, |y|<1, |2|<1 betrachtet, der zur Hälfte in Wasser von der Temperatur 0 
taucht und zur Hälfte an ein Medium mit linear ansteigender Temperatur grenzt. Es sei also 


Au=0 für |e|<1, |y|<1, 2]|<1 

u=0 für zs0 und |z|=1 oder |y|=1 

‘=2.für 220 und. Jel=1 oder RP N ee (18). 

u=0 für =—1 

ut fürzel 
Um nur einen ganz einfachen und dann allerdings auch noch recht ungenauen Ansatz zu beschrei- 
ben, werde 9, —=0, v,=1, 1y=2, 1 = x?-+ y?— 22? gesetzt, wobei die Symmetrien der Aufgabe 
benutzt sind. Die ganz kurze Rechnung liefert beim ersten Weg nach (12) die Gleichungen 


40 20 8:56 


2da,=8, Tr IANS a—=—4 


3 


und als Näherung 


1 1 15 
has + D5 2 118 (x + y—222) . 


1 
Am Rande schwankt v7 zwischen cr und eu und hier gilt im ganzen Bereich 
Ku 
336 


Beim zweiten Weg nach (17) hat man die Gleichungen (jetzt ist die Konstante zunächst 
fortzulassen) 


0,1172 Sup —us = 0,200. 


—8,=—4, —640,=8 
und somit 


1 1 ; 
„=0or72= 5 (@ + pP —22). 


°) A.G. Webster-G. Szegö, Partielle Differentialgleichungen der math tischen Physik, r 
und Bin. 1930, 8. 311. a N ematischen Physik, Lpz 
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Die Konstante ce kann man etwa so bestimmen: 9» —c— u schwankt zwischen den Werten 


1 5 3 3 

3 EI d AL . a ke . ER A . . 
EHER N RE setzen wir 6= 75,50 schwankt vo — u En 7 ‚und 16 gilt dann im 
e ganzen Bereich 3 FR 
4 I > \ [ur _ ul s 16 — 0,1875 0: ER eo = Pe (19). 

| x ® Die Fehlerschranken sind hier noch so groß wegen der Knickstellen in den Randwerten 
von u. Durch Hinzunahme weiterer Glieder im Ansatz lassen sich natürlich die Schranken leicht 


_ — —  herunterdrücken. 
Br: - II. An dem folgenden Beispiel sei gezeigt, daß man oft mit Hilfe einer Einpassungsmethode 
‚sehr rasch zur Aufstellung einer groben Näherungslösung gelangen kann. Es sei 


tn tut tu 0<esl, Wed... 


und u=s am ‚Rande = abe gehe und W—- Z) 


vorgelegt; läßt sich denten als Temperaturverteilung in einem Körper oder prismatischen Raum- 

_ teil von quadratischem Querschnitt in einem linearen Temperaturfeld u=x, während die innere 

- Wärmeleitfähigkeit linear mit & veränderlich ist. Durch spezielle Ansätze (u als Summe oder 
= a Produkt einer Funktion von # und einer Funktion von Y nn stellt man leicht spezielle Lö- 
Br „sungen % von (20) auf: 


ne 


e  a=l, y=h(l+a), weit, u=N, N (o+ 2) CRY, 


wobei Ns die Neumannsche A nullter aa: ist. Die noch freie Konstante ß 


a Bein etwas. re 
obi nsthae ungleich un ans 


FF BildT, 
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1 R 
Für den Quadratmittelpunkt 0-7 y=® z.B. erhält man so 


1 
0,5436 <ulz 0) < 0,5521 oder (>, 0) —= 0,5479 + 0,0043 . 


7 


6. Verfahren gleichzeitiger Annäherung von Differentialgleichung und Randbedingungen 
Es liege eine Randwertaufgabe für eine Funktion v'von n unabhängigen Veränderlichen 


%s...,%, vor. Die Differentialgleichung sei nichtlinear 
ou ou u Mu 
Ws end a EUER 21) 
a FI dar’ Da dee ) ( 
oder linear 
oArr ...+ &n h 
L[u]= Au. na rl U) DI (22), 
RR PR a le 


wo L[u]ein linearer homogener Differentialausdruck in u mit gegebenen, von den x; abhängigen 
Koeffizienten A, ... a, Ist. Es kommen nichtlineare Randbedingungen 


| ou ou u 0:u ) 
DR Re VE E 


a ton. dar om 
oder lineare Randbedingungen 


=0' mir Tr 
Kies) 


ou ou Hu ou 
De a ee er ae re 24 
.[ 5 0 0x. 022° 0x, 02, ) Pu et: Kt Sr 
hinzu, wobei U, gegebene lineare homogene Differentialausdrücke in v sind.‘ Dabei ist ® ein ge- » 
. gebener Bereich im 2,-2-... &,-Raum und 7’, sind (n—1)-dimensionale Hyperflächen; ge- 
wöhnlich werden es Randflächen von ® sein. Bisher wurden gewöhnlich zwei Arten von 


j 
Er Näherungsverfahren untersucht, bei denen « durch einen Näherungsausdruck w mit p Para- 
A metern a, ..., 


RR 1 RE ea RER Be SE nr ARD 
angenähert wird: 

1. Fall (Randmethode): w erfüllt für beliebige Werte a, exakt die Differentialgleichung 
(21); die a, sind so zu bestimmen, daß die Randbedingungen (23) möglichst gut befriedigt werden. 

2. Fall (Gebietsmethode): w genügt für beliebige Werte der a, den Randbedingungen 
(23); die a, sind so zu wählen, daß die Differentialgleichung (21) möglichst gut erfüllt wird. 

In beiden Fällen pflegte man die Güte der Näherungen nach der Größe des Fehlers &, d.h. 
der Abweichung w— u der Näherung w von der exakten Lösung « zu beurteilen und fragte nach 
Abschätzungen für den Betrag je| des Fehlers. Nun betont R. v. Mises?), daß es vom physi- 
kalischen Standpunkt aus vernünftiger ist, Toleranzen 6, ö, vorzuschreiben und nach einer Nähe- 
Tung w zu fragen, die diese Toleranzen bei der Differentialgleichung und den Randbedingungen 
einhält (in wohl leicht verständlicher Schreibweise) : 


IF[wl| <ö, A ir 


Bei physikalischen und technischen Aufgaben werden gewöhnlich weder die Differential- 
gleichung noch die Randbedingungen exakt gelten, sondern es sind in der mathematischen For- 
mulierung vereinfachende Annahmen, Linearisierungen und andere Idealisierungen enthalten, 
so daß (21) und (23) nicht streng gelten werden. 

Einen naheliegenden Ansatz für w erhält man aus der Forderung, den Integralausdruck 


k 
J= [Fol det Zku /, Vf lafss. „lie oe (27) 


zum Minimum zu machen. Dabei ist dr das Volumenelement in ® und df das Oberflächenelement 
auf J',, die Zahlen k, sind gewisse positive Konstanten, über deren geeignete Wahl in Nr. 9 ge- 
sprochen wird. Etwas allgemeiner könnte man 


k 
J*= I p(z;) F?[w] dc+ ku] ?,(%;)V&[w] d‘=Minimum . ...... (28) 
„1 Ty 
ansetzen, wobei p und p, als Gewichtsfunktionen positive Ortsfunktionen sind. 


dh iv R.v. Mises, On Network Methods in Conformal Mapping and in Related Problems (ent dem- 
cnst), . 
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Wir setzen hier p=p,=1. Die Forderungen für ein Minimum von J 


0J 


a =1..,P) s:. (29) 


dienen.dann zur Bestimmung der Konstanten a,,..., Ip: 

Setzt man die so ermittelten a, in die Näherung w und diese in (21), (23)-ein, so sieht man, 
welchen Fehler man noch bei F und V, hat und kann notfalls einen Ansatz mit weiteren Kon- 
stanten a, durchführen oder die Gewichte %, bei den Randbedingungen ändern. 

Hat man eine Näherung w gefunden, welche die Toleranzforderungen (26) erfüllt, so tritt 
die Frage nach der Fehlerabschätzung jetzt in der folgenden Form auf: Wie groß ist bei gegebenen 
ö, ö„ der Spielraum für u? D.h. man hat also eine Näherungslösung w und die gesuchte Funktion 
u, welche beide den Ungleichungen (26) genügen 


Renee lol... rw, (30) 
und man fragt nach einer Schranke S für den Betrag des möglichen Unterschiedes 
lel=|w—u|ss$ Ba FE ER Me (31), 


wobei S von den ö, ö, abhängt. 
Bei linearen Randwertaufgaben (22), (24) folgen aus 
Klul—r|sd, Lul—risd, U,Wwl—y|ls%, |U,lu]—ylsd, . - 82) 
für die mögliche Abweichung e die Ungleichungen 
PrENE 2 29., 7° HERE 20, ER (33). 
Hieraus ist die Schranke S für |e| in (31) festzulegen. 


7. Gleichungen des Näherungsverfahrens bei linearen Randwertaufgaben 
Der hier gegebene Ansatz ist auch für die Auffassung verwendbar, daß man etwa die Diffe- 
rentialgleichung (22) und die Randbedingungen (24) als streng und (27) und (29) als Vorschrift 
für ein Näherungsverfahren ansieht, welches sich von den sonst üblichen. darin unterscheidet, 
daß man sich beim Ansatz für w noch nicht um die Differentialgleichung oder die Randbedin- 
gungen zu kümmern braucht. Das kann rechnerisch ein Vorteil sein, weil die strenge Erfüllung der 
Differentialgleichung oder der Randbedingungen mitunter zu schwerfälligeren Ansätzen und 
damit zu mühsameren Rechnungen führt. In (33) fällt dann der Faktor 2 fort: 
Baledrr el d et NER (34). 
Macht man für w einen linearen Ansatz 


so folgen aus ; } N \ 

ein Ban,le BES au. ar 0006) 

Mm Me | si BI oE, o=1 
für die a, die linearen Bestimmungsgleichungen mit symmetrischer Koeffizientenmatrix: 

k 
JtElw]—r)Zim]drt 2 h SU Uteldf=0.... >. (37) 
oder Be | 
k k 
| J4twlEtw] de > ku [Ute] U, [usa = rL [w,] dr+ P3 kulYu U,[w]df (38). 
o=1 B u=1 u A= m 
(0=1,.759) 


8. Fehlerabsehätzung 


Zur Durchführung einer Abschätzung des Fehlerbetrages |e| aus (33) oder (34) — wir können 
uns etwa auf (34) beschränken — legen wir eine speziellere Aufgabenklasse zugrunde, und zwar 


die erste Randwertaufgabe bei der elliptischen Differentialgleichung (1) (2) 


! n 02u n% ou 
s LE ee Die ES OUZT. n Zihtl en (39). 
[ ] 22 ik 0%; 08% & 3 0%; 


"Dabei sind a,,, bj, c, r im Bereich ® gegebene Funktionen von 2, .. ., %- Die als symmetrisch 


angenommene Matrix der a,, sei in ®-+-I" positiv definit, und es sei c6>0. Auf dem Rande 7" 


von ® seien für u Randwerte vorgeschrieben. Es ist also der Fehler e abzuschätzen aus 


eg Beau I nee (40). 


nn ist der Satz vom Randmaximum aus Nr. 2 anwendbar. 
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Nun sei 2 eine Funktion mit r 
Lej<—4<0n BREI (41). 


Oft wird man eine Funktion z mit L[2] = — 4A <0 sehr leicht angeben können, vgl.Nr.3, z.B. 


0 A 4 ET 


durch geschickte Wahl von z läßt sich die Fehlerabschätzung verbessern. Auf dem Rande 7’ 
sei 2, <2<i2,. Die Zahlen z, und z, sind als bekannt anzusehen. Zunächst sei c=0. 
Dann erfüllt nach (34) die Funktion 


9=:+ ut. Digi Bee (42) 
die Voraussetzungen des Satzes vom Randmaximum, 
Auf dem Rande ist 9 <sö,+ en also gilt diese Abschätzung auch in ganz ® und damit 
3 N 
Zr) a RE (43). 
Entsprechend ist Z[9] Z 0 für die Funktion 9= e— 7 ö; das Minimum von 9 wird auf /’ an- 
genommen, also ist 
DE OS. 
P 1 “7 2 
Somit gilt in Bez — dı— 3 (2,—z) und mit (43) insgesamt 
E; 0) 
TR lel sdı+g (BB) 5 u ee a ER (44). 
Im Falle c> 0 wird das Maximum von |p] am Rande angenommen, also nach (42), wenn 
| auf 7’ gilt 4 <B, 
Rn: le+22 <ö,+öB. 
se und damit E 
Er Fr le| = 6446 (B+ MM ) sö,+2 öB Ten a ee ie We (45). 
y # * u 
Fi 9. Beispiel 
B- Es sei 


—Au=0in B, u(z, +1)=0 für |e| sl, +1, y)=1—y für |y|sıi . (46) 


vorgelegt; Bist dabei das Quadrat |e|<1, |y|< 1, u(z, y) läßt sich deuten als stationäre Tem- 
peraturverteilung in einer quadratischen Platte mit gegebenen Randtemperaturen. 


In dem Ansatz, in dem Q das Quadrat O<2<1, 0 <y<{1 bedeutet und die Symmetrien 
der Aufgabe berücksichtigt sind, 


1 1 
I=[ftawrazaytk] [w(z, 1)]?de + H [w(1, y+y—1]? ay) ARTE TAT 
ist zunächst das Gewicht % geeignet zu wählen. Eine Hilfsfunktion 2 nach (41) ist hier z.B. 
- (2?+y) mit 4=1. 


1 1 ? 
Wegen c=0, n=5» 0<s %—2Sz lautet (44) Je] <s4+50, d.h. ein Fehler ö in der Dif- 


ferentialgleichung geht hier mit dem halben Gewicht wie ein Fehler ö,in der Randbedingung in 
die Fehlerschranke |e| ein; wir wählen daher k—=2. 
Für w werde ein Ansatz gemacht 


va, +0; Y’+a, (a! + y)-+a, 28 yR. 


(Es würde vielleicht näherliegen, a,2°-+a,y* an Stelle von a,(x*+y) zu schreiben, aber es er- 
gibt sich sogleich a,—a,); man kann nun w in (47) einsetzen und die Gleichungen (29) bilden; 
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etwas bequemer ist es, direkt die Gleichungen (37) aufzustellen: 


0 20,1 2004 360,1 10,—= 10 
140a,+ 3360, 280a4;,+1016a,+ 1960,— 70 
140@,+ 2800,+ 3360,+10160,+ 1960,— 14 
378a,+1620a,-+1428a,-+ 75324, 13260, — 210 
2100, + 5880,+ 5880,+2652a,+ 5180;— 42 
mit der Lösung 
a1= u (0,0)= 
a=0,50011, a,——0,49989, a,=0,096%, a,—=—0,58196. 
Für die so bestimmte Funktion w wird 


Aw=2(a,+a;)+2(60,+4;) (#24 y?)= en 2 


— 0,41086, 


es ist daher —0,00121 < Aw< 0,00044 


j 22 
i dw| <ö 1gog7 > 9001203 . 
‘ Ferner brauchen wir zur Fehlerabschätzung den Randfehler; es ist 


1 
vw, 1)=, ++, +(+ a,)0?+a, 2! = 3. 18981 (1155—11970 «+ 14175 x*) 


vll, y)+y’—1= En (1155 — 11970 y? + 14175 9). 


420 


Der Fehler e=w— u liegt am Rände zwischen den Schranken —0,0094 und 18981 


für die Funktion w* = w-+-0,0068 würde e&®=w*— uam Rande zwischen 
—0,0162 und —-0,0162 

schwanken, so daß man ö, = 0,0162 und damit für |e*| die für den ganzen Bereich gültige Schranke 

0,0162-4 = 0,00121 oder 0,0169 (man könnte auch 0,0168 erreichen) erhielte. Für vergeben sich 


damit die Schranken 


—=(0,02298; 


w— 0,0237 su sw+0,0101. 


Diese Fehlerabschätzung gilt für das ganze Quadrat und gibt die Größenordnung des Fehlers 
richtig wieder, da ja in der Ecke tatsächlich w um 0,02298 von u abweicht. (In den Ecken liegt 
eine Singularität von u vor, da wegen .der Randwerte dort 


Va 
An 2) 


ist, während im Felde Au=0 sein soll). 


: Eingegangen am 1. Oktober 1951. 


Über die Entstehung von Ringwellen an einer Flüssigkeits- 
oberfläche durch unter dieser gelegene, kugelige 
periodische. Quellensysteme*. 

Von Hsien Chih Liu in Schantung, China 


In einer idealen, nur der Wirkung der Erdschwere unterworfenen Flüssigkeit mit unendlich ausgedehnter 
freier Oberfläche und unbegrenzter Tiefe wird eine unter der a N oBerfläche gelegene, kugelige perio- 
dische Quelle als Störungszentrum angenommen. Für den Fall, daß die Zeit über allen Grenzen wächst, 
wird die dadurch hervorgerufene Oberflächenerhebung aus dem Geschwindigkeitspotential der Quelle auf- 
gestellt. Das Verfahren wurde auch auf den Fall angewandt, daß die Ergiebigkeit der Quelle über einer senk- 
rechten Strecke gleichförmig verteilt ist. 

In an ideal liquid with infinitely extended free surface and infinite depth, a spherical source, situated 
below the surface, is acting periodically as centre of perturbation. The only external force is gravity. The 
resulting displacement of ihe surface is computed from the velocity potential of the source, öf the time co- 
ordinate is indefinitely increased. The method was also applied to the case, that ihe source is equally extended 
to a straight line, normal to the surface. _ = 


*) Aus einer von der Mathematik-Naturwissenschaftlichen Fakultät der Georg-August-Universität 
zu Göttingen angenommenen Dissertation. 
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Dans un liquide imaginaire, soumis seulement & l’influence de la pesanteur de la terre, ü surface libre et 
stendue & l’infini et & profondeur illimitee, une source periodique de la forme d’une boule situee au-dessous de 
la surface du liquide est suppos&e comme centre de trouble. Au cas que le temps croßt au-delä de toutes limites, 
l’&lvation de la surface causse par cet accroissement, est constitude du potentiel de velocite de la source. La 
methode fut aussi appliqude au cas que le rendement de la source est distribus &galement sur une &tendue 
verticale. 

B uneanbHof 3KHNKOCTH, HaXonammelich HOA BIHAHHEeM TOJIBKA 3EMHOTO IPHTAKEHHA, C HEOTPA- 
HHYeHHOl CBOÖONHOH IOBEePXHOCTBO U HeOrPAHHYeHHOH T.LYOHPE, IPHHHMAETCH B BIaYeCTBe 
BO3MYINAlINeLO NEHTPa chepmyeckuä, ePHONHYeCKHÄ UCTOYHHK, P&CHON0KeHHLBIÄ ION HOBePX- 
HOCTbIO 3KHAKOCTH. BblsbIBAeMOe 3THM BO3BEIIMEHHe TIOBEPXHOCTH YCTAH&BJINBAETCH H3 TIOTEH- 
IHaa CKOPOCTH UCTOYHHRA Ha cıIyyali ÖecmpemeiıbHOrO BO3PACTAHUA BPEMEHH. Cnocoß npy- 
MEHALICH H K TOMY CIIY4A, KOFAa SbheKTHBHOCTL UCTOYHUKA PABHOMePHO PäacıpepeleHa Ha 
BEPTHKAJIBHbEIM OTPE3KOM. 


1. Einleitung 


Mit der Entstehung von Ringwellen durch einen kurzen Anstoß auf eine Flüssigkeitsober- 
fläche haben sich zuerst Cauchy!) und Poisson befaßt. Dabei wurde die Störungsstelle als 
unendlich klein angenommen; die Stoßkraft ist in ihrer Größe jedoch beschränkt. Wegen der 
großen mathematischen Schwierigkeit wurde die Behandlung noch weiter spezialisiert. Es wurde 
vorausgesetzt, daß einerseits die. Wirkungsdauer der Stoßkraft als sehr klein betrachtet wird 
und andererseits die Wirkungsdauer der Stoßkraft als sehr groß zugelassen ist, aber nur die 
Wellen in großer Entfernung von der Störungsstelle studiert werden. Dieselbe Untersuchung 
haben später Widenbauer?) und Hinze?) hinsichtlich der mathematischen Behandlung 
mit geringeren Einschränkungen der Bedingungen wesentlich erweitert. 

Über die Erzeugung von Oberflächenwellen durch eine unter der Flüssigkeitsoberfläche 
befindliche Störung hat Holstein) sich mit einer horizontal liegenden pulsierenden Quell- 
linie beschäftigt. In seiner Ausführung wurde die Behandlung auf den zweidimensionalen Fall 
beschränkt. Trotz der geraden Front der so entstehenden Wellenzüge hat seine Untersuchung 
hinsichtlich der Art der Entstehungsursache doch enge Verwandtschaft mit der vorliegenden 
Arbeit. Für einen Vergleich der Wirkung von Störungen auf und unter der Flüssigkeits- 
oberfläche seien noch die Veröffentlichungen von Schuler) und Dimpker®) erwähnt. 


2. Voraussetzungen, Randbedingungen und Problemstellung 


Im Falle einer wirbelfreien, nur der Wirkung konservativer Kräfte unterworfenen Strömung 
geht man von der Bernoullischen Gleichung aus: 


A: vo  0® 
\ statt 


Dabei bedeutet: ”—=v,-+vy+v; — Quadrat des Geschwindigkeitsbetrages, ® das Geschwindig- 
keitspotential, U das Kraftpotential, p der Druck, o die Dichte und F(t) eine willkürliche Funktion 
der Zeit. Wenn man die willkürliche Funktion F(t) in dem Wert von d®/öt einbegriffen be- 
trachtet und sich auf kleine Ausschläge der Flüssigkeitsteilchen beschränkt, so daß v? in Glei-. 
chung (1) vernachlässigt werden kann, folgt weiterhin mit U=—gz 


er od 
2 ; Dabei bezeichnet z die Koordinate, die in die Richtung der Erdbeschleunigung gfällt. Bezeichnet 
Bi man nun mit £ die Erhebung der Oberfläche über den ungestörten Spiegel, so folgt aus (2) mit 
2 Y z={ und p=0 (da der atmosphärische Luftdruck an der Oberfläche überall derselbe ist) 

Be VE Ne | 
. = Far] (gültig Tür, ==), 2.10 Sl Bee (3), 
n Bi die angenähert auch für 2=0 gilt. 

Ber. ' Weiter darf bei den vorausgesetzten kleinen Amplitud i Krümn 
3 A plituden bzw. kleinen Krümmungen der - 
ER %& 98 F 
67 3 (gültig fürz=0)....... (4) 
Bi EN 

en ) H. Lamb, Lehrbuch der H d. „We Fec 
4 2) Z.angew. Math. Mech. 14 ae, Kal. ee 


®) Z.angew. Math. Mech. 16 (1936), 8.143. . 
4) Z.angew. Math. Mech. 17 (1937), 8.38. 
5) Z.angew. Math. Mech. 16 (1936), $. 65. 
®) Annalen der Physik 19 (1934), S. 225. 


“ die durch eine, solche periodische Quelle hervorgerufene Ober- 


Oberfläche besitzt und der Wirkung der Erdschwere unterliegt. 
Die Flüssigkeitsbewegung in unmittelbarer Umgebung der 


im übrigen Flüssigkeitsbereich noch zusätzliche Potentiale 
$ ‚gendermaßen. fest (Bild 1): Die r-Achse liegt in der freien 


‚außen. Die z-Achse weist in Richtung der Erdbeschleunigung g. 
_ Der Quellpunkt liegt bei r=0, z=h, wobei der Wert h kurz 


 gleichung die folgende Gestalt an 
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gesetzt werden, als Bedingung dafür, daß die Normalkomponente der Oberflächengeschwindigkeit' 


gleich der Normalkomponente der Teilchengeschwindigkeit an dieser Stelle sein muß. Führt 
man & aus (3) in (4) ein, so wird 
od 1098 
%2 gg 
Gleichung (5) ist bekannt als die Oberflächenbedingung, wenn die Oberflächenwellen im Vergleich 
zu ihrer Wellenlänge kleine Amplituden haben. 

Aus (3) ist ersichtlich, daß ein bestimmtes Geschwindigkeitspotential ® eine Oberflächen- 
erhebung hervorrufen kann. Daher besteht die Aufgabe dieser Arbeit einfach in der Bestim- 
mung solcher möglicher Geschwindigkeitspotentiale und damit auch des dadurch hervorgerufenen 
Wellenzugs für den dreidimensionalen Fall unter einem angegebenen Geschwindigkeitspotential 
als Ursache. 

Zur Vereinfachung des Problems wollen wir unsere Behandlung grundsätzlich auf die 
eingeschwungenen Bewegungen beschränken, Mathematisch ausgedrückt, müssen wir dann die 
endgültige Lösung unseres Potentialproblems in der Gestalt 

sin a 


N N 5 
= 3 dla, y,2,0)= 3 ple, y, 2) | (6) 
A=1 i=1 : cos wit 


ansetzen. Wir können verschiedene Potentiale der Störung Ö, erhalten für eine Störfre- 


(BBlllO 2in 2,0). 2 meld), 


 quenz . Diese muß außer der Oberflächenbedingung (5) noch verschiedenen Bedingungen 


genügen”). Die L’aplacesche Differentialgleichung = 
x DD OEDE OD. 7 
tar FE ET a Bw) 


muß im ganzen in Betracht kommenden Flüssigkeitsbereich erfüllt sein. Außerdem muß man 
noch fordern, daß die Flüssigkeitsbewegungen in großen Tiefen unbeschränkt abklingen. Das 
ergibt für die Vertikalgeschwindigkeit die Bedingung 


3. Eine kugelige, periodische Quelle als Entstehungsursache 
a)Die Bestimmung der Geschwindigkeitspotentiale. In einem 
allseitig unendlich ausgedehnten Flüssigkeitsbereich ist das Geschwindigkeitspotential einer 
kugeligen, periodischen Quelle gegeben durch: 


Ask, 


Punkt des Potentials ®,. Wir stellen uns nun die Aufgabe, 


flächenerhebung zu bestimmen, wenn die Flüssigkeit eine freie 


Quelle wird durch die Wirkung von ®, bestimmt sein, während 


existieren werden, deren Bestimmung unsere nächste Aufgabe 
sein wird. Wir legen ein zylindrisches Koordinatensystem fol- 


Oberfläche und weist von dem Nullpunkt r=0 radial nach 


mit „Eintauchtiefe der Quelle‘ benannt und als Parameter 


Bild. Festlegung des Koordinaten- 


betrachtet wird. Für den vom Quellpunkt ausgehenden Ra- systems (r,z) und der Quelle Qu im 
 diusvektor R, gilt offenbar die geometrische Beziehung Eieesa elaberaich 
en . Er I Er j Fa R} er Dar (z—h)? ee‘ Sr (10). 


Be Durch die Einführung eines zylindrischen Koordinatensystems nimmt die Kontinuitäts- 


ARE ne — ö ad = 
FE Eee > . Leer ee: (De ee (11). 
Er ara Mer? | 


15 


= —0 (gültig Tür’ 2.00, 20.70... (Oz 


COS RZ Be a Se Yo ee (I). 
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Die übrigen Gleichungen haben ihre ursprüngliche Gestalt beibehalten. Jedoch muß man jetzt 
statt ®,(z,y,2;t) und &(2,4; t) schreiben &;(r,2;t) und &(r;t). Zu 

Um die Randbedingungen an der Oberfläche zu erfüllen spiegeln wir die Quelle an der 
Oberfläche. Diese Quelle ist dann gegeben durch das Potential 


d,=— 


wobei BR Er Hl Er RT (13) 
gesetzt ist. Durch Einführung des zweiten Potentials , ist also die freie Flüssigkeitsoberfläche 
z=0 durch die folgenden Beziehungen ausgezeichnet 

8, +9,=0 (gültig für z=0) . - - . (14) 


und 
00, _ 0 ___ Ch _ 
%2 02 An(r?+A2)%' 
Mit Rücksicht auf (11) wählen wir als drittes Potential die unter dem Ausdruck ‚‚Potential 
der freien Welle‘ bekannte Form ®) 


2, = AJ (de: sin wre er ee (16), 


wobei A ein willkürlicher Faktor und J,(Ar) die Besselsche Funktion nullter Ordnung bedeuten. 
Das Einsetzen von ®, in (5) ergibt für A die Bestimmungsgleichung 


cos wt (gültig für z=0). . (15). 


damit ©, allen vorgeschriebenen Bedingungen genügt. 
Da die Oberflächenbedingung (5) noch nicht durch die Summe ®,—+-6, erfüllt ist, wählen. 
wir zweckmäßig ein weiteres Potential vom Typ ®,, aber nach Fourier verallgemeinert 


Be Ir Ju ler)» e FEdE: cos at Kuna Bean (18), 


worin die Funktion f(&) uns die erforderliche Anpassungsmöglichkeit bietet. Das Einsetzen der 
Summe ®,+®,-+-9, in die Oberflächenbedingung (5) liefert 


© 2 Q ; h 
| ( | FE)Is(lE r)dg gs PETCREWNT alle 05.2 2807 vepkel Fr Biker (19). 
e . 
Wenn man die bekannte Beziehung °) 
T 1 
J er J (Er)de = re] 


beiderseits nach dem Parameter h differenziert, was nach dem Leibnizschen Satz !°) durchaus 
erlaubt ist, dann erhält man 


ü h | 
[® En] (Er)&de = Thames (20). 
{ . 
Die Kombination der Gleichungen (19) und (20) liefert uns die‘gesuchte Funktion 
RT 
Se Ve =: (21) 
ESEF 
Mit Rücksicht auf Gleichung (18) wird schließlich das vierte Potential 
PER REE e | 
=.) ® PR) ler) da ee (22). 
a 


Jetzt sind alle Bedingungen durch die vier Potentiale erfüllt und damit kann das Potential- 


problem im Prinzip als gelöst betrachtet werden. Es bleibt i i 
ee ibt noch die Ausführung der Integration 


' *) H.Lamb, Lehrbuch der Hydrodynamik, 8.483. Teubner 1931. 


°) G.N. Watson: A treatise on the theory of Bessel functions. Cambri | 
ae, WR, Osgood: Lehrbuch der Funktionentheorie 1, 8.120. Teubner ee DaB IE: 
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b) Das Potential &, im Gebiet großer Entfernung. In Gleichung (22) wollen 
wir benennen 


oo 


Te ei SolcH de, a we ann 23 


Res 
a 


0 


Da wir zunächst das Potential ®, im Gebiet großer Entfernung vom Koordinatenursprung 
errechnen wollen, sind wir berechtigt, J,(&r) durch seine asymptotische Näherung*) zu ersetzen, 


somit gilt 
3 (lei 
=> e-h+2)8 a ea ER ee 2A), 
0 


Hierin führen wir statt & die neue Variable u—yVE ein und schreiben v„=h-+z— ır. Außer- 
dem setzen wir I gleich dem reellen Teil einer komplexen Funktion und somit wird 


SW SE 2 
1= Rel2 ei u De Er REN Br Bee a8 
) TTr 


Nun schreiben wir 


Denken wir uns nun in der komplexen u-Ebene eine lückenlose Kontur (Bild 2) zusammen- 


2 
gesetzt aus der reellen Achse von Null bis R, wobei R?> = ist und der singuläre Punkt u= = 


Y9 


mittels eines halben nach oben gerichteten Kreisbogens umgegangen wird, ferner aus einem 


.T ‚N 
A U ‘7 : : Be 
Kreisbogen von k bis Re * und einer Geraden, die von Re * bis 


Null zurückführt!!). Innerhalb des von dieser Kontur einge- Reit 
schlossenen Bereiches ist der Integrand regulär, und nach 
dem Gauchyschen Integralsatz ist das Integral an der Kon- 
tur entlang gleich Null. Auf dem Kreisbogen von R bis Re: 
bleibt Re(v u) >0, so daß der Betrag des zugehörigen Integrals 
mit Ro nach Null geht. Wie bekommen also für das 
Integral 
‚A 0 w 
©» 5 4 Ye R 
er, & u2 Bild 2. Integrationsweg in der 
I, —_ınT e En erw — 2 LEETERE (27) komplexen u-Ebene 
g er! 


In (27) stellt das erste Glied rechts das halbe Residuum für den Punkt w= © dar. In dem 


- 9 
letzten Integral von (27) können wir den Integrationsweg durch eine beliebige unter dem Winkel o 


von Null aus nach unendlich gezogene Gerade ersetzen, solange nur 0 <P<Z bleibt. Das Er- 


BaRo= = w® 
gebnis wird unabhängig von 9 sein. Schreiben wir deshalb “= = e 2. .gy=o, 7 —=ß? und 


* Absehhtkumn des Restgliedes wird im Abschnitt (e) ausgeführt. 

2) Vgl. Frank und Minen. Die Differential- und Integralgleichungen der Mechanik und Physik, 
1, 2. Aufl. (1930), 8. 145. 
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setzt man weiter &ß?=4#?, so bekommt man 


. @ er 2 
I, =ın — e 0 +.— ee" 


2yY Ya 


0 
Die Integration kann ausgeführt werden und liefert: 


= a he ee In j2n+2 . (29). 
heine" (in ra Ian a | 


Wenn », « und Ain I, beseitigt werden, erhält man aus (25), (26) und (29) 
I= Re: bio] Pe 7 
r 8 
RER 1 31 
2 (42m) ) 


1 nn I TI FE 
3” erde *.335- Qu) 


Kehren wir zu unseren Ausgangsbezeichnungen zurück, so erhalten wir für das Potential &, 
in großer Entfernung von dem Koordinatenursprung schließlich: 


4 t BR NEN Er 

® Q 2 Vz le cos ur: A, RIES TREUL IE. "hte 3 "htz 

= — 120 |)/ —: —r— — 
2nfr 9 g g [(R+2)? +} 


(2 var)” 


1 1 
cos. arctg a sin — erer 


au [(h+2)?-Fr?% Die Su an .(2n-+1) (31). 
: 1 r 
f sinzarctg 7, — COszareig7 
x cos (Bl) area er ae 


2 n+1 
> (2 ar) 5 
—1)r+! — I —— ——— sin (n+1) arcetg-——[ cos wt 
EN ) 1-3-5---(&nt1) RT NARTTT 
Hier bedeutet der Index g von D®,, das Gebiet großer Entfernung von dem Koordinatenursprung. 
c) Die stehenden Wellen im Gebiet großer Entfernung. Wenn (31) in (3) 
eingesetzt und die Bedingung z=0 beachtet wird, dann ergibt die Gleichung der stehenden 
Wellen in großem Abstand von dem Koordinatenursprung 


1 1 
Eat Be sin -arctg— „+ 008; aretg 
u=-— a v- naar“ n.e Band ( r+- 2)+ Se 
rt ; +] 
ae N N 
sin —aretg” Hooszareigt „ 2ayı+@)) 
> a RR RR ER BT 
1-35. -(2r+1) 


Il 
o 


BER 


RS AU 68 ß ATCT E 
SEIET 3 DT, 


SE 
Say 27) 


{ r\., 
13-5. (2n +1) -sin(R-+1) arctg—|sinwt. 


(32). 


x cos (n+ 1) arctg z + 


3 
ll 
o 


ee nenn ee En SW Zu ze ws ERSTE TEUER 


ee I Br 'W 
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Es ist leicht nachzuprüfen, daß ®, und Ö, an der Oberflächenerhebung der Flüssigkeit 
nicht direkt beteiligt sind. 


Für das Gebiet von sehr großer Entfernung kann gegenüber r vernachlässigt werden, 
was selbstverständlich ist, weil für die Entstehung bemerkbarer Wellenausschläge der Wert 


der Eintauchtiefe } nach oben begrenzt sein muß. Damit erhält man mit arctg > zde 
folgende einfachere Formel 


Qw | ae er =) 
(En 2\azre Seas; | rer) 
. zgrV2 | 9) y & 4 u 
2 \n+1 

: 07 | 
— —1)®+1 : cos (n+ 1) —(sin wi 

> Teer 3] 
Jetzt kehren wir zur Diskussion der Gleichung (32) zurück. Gleichung (32) besteht wesent- 
lich aus zwei Multiplikationsfaktoren, von denen der eine nur von der räumlichen Veränder- 


lichen r, der zweite nur von der Zeit abhängt, und zwar so, daß das Produkt zeitlich periodisch 
ist und bei der Deutung von Z,, als Verrückung harmonische Bewegungen der einzelnen Punkte 


83), 


‘ des in Frage kommenden Raumbereichs darstellt. Zur graphischen Darstellung von (32) ist es 


praktisch, sich vorzustellen, daß die kugelige Quelle um eine mittlere Amplitude mit der Kreis- 
frequenz ® atmet. Um ein Momentbild herzustellen, setzen wir willkürlich der Beseitigung 


des negativen Vorzeichens zuliebe ai . Der Bequemlichkeit halber zerlegen wir (32) in 
die folgenden einzelnen Ausdrücke 
Dee w* 4 
el ) 
uno 


7 AngVhr 


ka" 


| da 1 r 
sin 5 arcig ;R —+ cos — arctg ya 


— 2 S n+1 3 7 
x SWESPFEFTNDDIITEE 122041) = nel) cos (n-+-1)aretg 3 
BB) 2 
DEN 
RL 
= 2ngVhr 
2 seo 
sin > areig 4 — cos Zaretgt 2 ( 5 "Vı4(g) ı 
ee Da er 1 Eure a a 
| 146) = | 


Für unsere numerische Auswertung setzen wir willkürlich die Ergiebigkeit der Quelle 
Q=10 Liter/sek, die Binfeuchtiefe h=10 cm und ae 1. Als Bereich sei noch 1 bis 11 
zu wählen. Die errechneten Ergebnisse sind in Bild 3 zusammen graphisch dargestellt. 

Der Faktor cos (+2) in n, liefert uns eine einfachharmonische Schwingung. Die 


Knotenpunkte bestimmt man aus der Beziehung 
ZEREER > a ON ent (34) 
w? 4 


und die Wellenlänge beträgt 


Ihre Länge ist umgekehrt mit dem Quadrat der Kreisfrequenz & proportional. Für unser Bei- 


spiel ist A=62,80 cm. Das Erscheinen von Yr im Nenner bedingt eine mit Zuwachs der Ent- 
fernung r abnehmende Amplitude für n,. Wie die Kurve in Bild 3 zeigt, liefert n, keine einfach- 


« 


x 2. . Math. Mech. 
218 Hsien Chih Liu, Über die Entstehung von Ringwellen 34.32 Xr.7 Juli 1958 
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1 j j i ößere Länge als die 
harmonische Schwingung und zwar hat die untere Halbschwingung eine größ: 
obere, dagegen sieht die obere Halbwelle spitzer. aus. Die Amplitude klingt mit Zunahme von r 
auch ab. Das dritte Glied n, ist gegenüber n, und 7, vernachlässigbar klein, infolgedessen ist 


Bild 3. Momentbilder für die einzelnen Oberflächenerhebungen n,, ns und n,in Abhängigkeit von der Entfernungr 
für den Zeitpunkt wt= 3 n 


das Verschwinden von n; in (33) bei sehr großer Entfernung tatsächlich eine richtige Folgerung. 
Für Abnahme der Amplitude spielt r hier dieselbe Rolle wie bei n, und n,. Die obere und untere 
Halbwelle von n; haben auch ungleiche Länge wie bei n3. 


Die graphische Darstellung von &,,=n1ı+Na+n3 findet man in Bild 4 vor allem mit einem 
anderen Maßstab für die Ordinate. Obgleich die obere und untere Halbwelle keinen großen 
Unterschied in ihrer Länge haben, ist die Schwingung keine einfachharmonische, da die Kurve 


noch durch das Auftreten von /r im Nenner der Wellengleichung verzerrt ist. 


NN) 
Kar 4 
TERN BZ EEE 
REG 


Bild 4. Momentbild für die Oberflächenerhebungen 7, ,n, und sin Abhängigkeit von der Entfernung r für den Zeitpunkt wt= Er 


Die errechneten Ergebnisse bekräftigen ferner, daß n, in &,, seiner Größe nach erheblich _ 
dominiert. ‚Bei vorgegebener Ergiebigkeit und Eintauchtiefe der Quelle ist es von starkem 
Interesse, eine passende Kreisfrequenz so festzustellen, daß n, seinen Größtwert annimmt. Es 


: A 


. drat der Kreisfrequenz der Quelle wird. 
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ist leicht festzustellen, daß für gegebene Eintauchtiefe A 


Sa ee 


die Bedingung für ® ist, bei der 7, sein Maximum findet, und zwar wird der maßgebende Faktor 
von Nımas nach vorstehender Gleichung für n, 


SUN ) 
(are : N eh‘ 


Ferner sei noch darauf hingewiesen, daß, wenn der Quellpunkt an die freie Oberfläche rückt, 
die Wellenamplitude von », in großer Entfernung bis auf den Faktor FE proportional dem Qua- 
r 

d) Die fortschreitenden Wellen im Gebiet großer Entfernung. Die Rechnung 
im Abschnitt c) hat gezeigt, daß n, und n, mit Zuwachs von r rasch abklingen. Daher herrscht 


im Gebiet sehr großer Entfernung von dem Koordinatenursprung tatsächlich nur die einfachere 
stehende Welle 


2 we, 20% w? 7% 
GA ee ee g . Bene .<1 
Zum —Q wor g e cos (® +2) SID. 0. Ser la: 


2 
Wenn man 9, aus (16) in (2) einsetzt und wiederum J, e ) durch seine asymptotische 


Näherung ersetzt, dann hat man . 


ER 
wobei der willkürliche Faktor A in ®, durch Dimensionsbetrachtung und die notwendige, 
exponentielle Wirkung der Eintauchtiefe 3 bestimmt ist. Nach Kombination von (37) und (38) 
erhält man : 
Die w@® 7 
— N ER ET Be Er ART : 
Gmb t Lu Q na cos (*r 1 0) (39) 
Dieser Ausdruck repräsentiert einen fortlaufenden Wellenzug. Die Amplitude der Welle nimmt 
mit Zuwachs von r stets ab. 


e) Abschätzung des Restgliedes von ®,für große Entfernung. Im Abschnitt b) 


wurde die Besselsche Funktion J,(&r) durch ihre asymptotische Näherung Vz- :coS \: —i) 


ersetzt. Da diese Ersetzung unter dem Integral vorgenommen wird, und sich das Integrations- 
gebiet von Null an erstreckt, muß noch gezeigt werden, daß diese asymptotische Annäherung 
im vorliegenden Falle erlaubt ist. 

Um die Abschätzung bequem ausführen zu können, macht man zweckmäßig die Integra- 
tionsvariable durch die Substitution £r=e in dem Ausdruck 


oo 


I= | etta: = -Io(Er) de 


Pe 


9 


dimensionslos und erhält somit 


ale. 


Für großes Argument e kann man schreiben 


ee - Vz Ei (3) +0(le}>®), 


und somit bedeutet die Restabschätzung, die wir gleich ausführen werden, einfach die Fest- 
stellung, ob mit Zuwachs von r und für das IntervalO se<w eine genügend kleine obere 
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Schranke für den Wert 


ee £) 
142, Jo(e) — V- - cos (e-2) 
r-Al= et —— de 
er 
0 ge 


existiert oder nicht. Anders gesagt muß man verlangen, daß auch rAI mit r—> co gegen Null 
geht, da I selbst mit 1/r gegen Null geht. 
Den Integrationsweg wollen wir in der komplexen 


Sr e-Ebene (Bild 5) so wählen, daß er längs der reellen Achse 
Pr © von Null 5 Mr nuBe erstreckt wird, wobei der singuläre 

-/@--- 
0 2 Sr 2/ Yr- 1) Punkt e= = r mittels eines halben in dem ersten Quadran- 


Bild 5. Integrationaweg in der komplexen. ten liegenden Kreises umgegangen wird, wie die gestrichelte 
Das Linie zeigt. Für unsere Abschätzung ist es zweckmäßig das 
Integrationsintervall so zu zerlegen, daß er aus den Anteilen 


2 2 
VER Br2rEr e ‚—1) [über den Halbkreis]und 2 Ü —) ... 00 besteht. Diese Zerlegung 


bewirkt hauptsächlich, daß in dem ersten Intervall O...2 sowohl J,(e) als auch cos a 
positiv bleiben. Somit erhalten wir jetzt ein ae für rAI 


rn) 
r- AI=r(AlI,A +AlL,+Al,)) = f I» EN: ee odle NOReR? 
En 


= (ri) ge 
2 
| Außerdem muß noch gezeigt werden, daß r-AI im Pol 3: auch beschränkt ist. 
E. &) Abschätzung des ersten Anteils rA/,. Da der Wert : — in dem Inter- 
SE w“r 
BR ee 
= vall Ose<s2 sich monoton ändert und 
1 €g (eg \* 
ar Ef 
a Pk) 
ge 
ist, gilt 
1 29  /2g\" 
ui) 
ge 


Da ferner in diesem Intervall sowohl die Besselsche Funktion J,(e) als auch cos (3) 
positiv sind, so folgt offenbar die Ungleichung 


5 EaE Ba 
20-2, cos (e-3) </z cos (12)< = ; 


Daher kann man schreiben 


2 — 
BR 7 2 
n+z  Jo(E) N — (cos (e-2) 2 
a Bere TE 4 3 Ya“ ht. 
r AL|= BALE Be ARE NEIR Tee 2: <a 1 2 (28) = [- r ie 


wer in LAT, 7 
0 : ge 0 & 

2 

de 4 Tre 

Mit ra 2Y2 und —— 57’ wobei A die Wellenlänge bedeutet, erhält man 
4 4 
|AT, | a E A ’ 
nr ei 


Z. angew. Math. Mech. ‚ . . . n : 
Bd.32 Nr.7 Jui1ıg2 Hsien Chih Liu, Über dis Entstehung von Ringwellen 221 


wobei die kleine Größe höherer Ordnung, die unter dem Summenzeichen steht, weggelassen 
wurde. Also ist r AI, nach oben beschränkt und nimmt rasch ab, wenn r>A ist. 


ß) Abschätzung des zweiten Anteils rA],. Man kann leicht zeigen, daß die Un- 
gleichung gilt 


(771) = 7 
hte, J,(E) — ee ee 
RER = r de 


e 


ee 
2 gr 
2 (© r-) 
ent), 
= air 02,00 e—7))de. 
2 g2 


Da der Ausdruck in der Klammer längs des Halbkreises beschränkt ist, so folgt 
Is 
1 23 ( z 
|rAT,| ST of [0-] BE Er SH )ae. 
ber 


mn 


Demnach geht rAI, mit r— co nach Null. 
2 
y) Abschätzung des dritten Anteils rAI,. Längs der reellen Achse für >2(“ -i) 


und r> I wird 


@® 
-_ = 
ern 
ge 
Daher hat man 
? h-+z FEN 
— 2 PyA 
rAT,| z2 j (RAR (ne Ve (Zar. 


er 
Wendet man wieder die asymptotische Abschätzung 
DE 7 
J, (e) Vz cos (e= a) < A(le|?) 


an, welche für hinreichend große e gültig ist, so hat man 


oo 


h+z de 


var =24 | ae - 


en) 


f) 
Wenn man a e=n Setzt, so ergibt sich 
r 

h+2 ! dn 
AI 24 ee : 

a FE nn 

wir = h42 
2(, r- ) n 


1 
Demnach geht rA/; mit — gegen Null. 
Ir 


Ir 
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RE ae! nl is es 


2 ; 
6) Abschätzung des Restglieds bei dem Pol 7 Das Halbresiduum wird 


an w®? 29 w*? 7 
Te (6 Var 5): 


welches gemäß der Restabschätzung für die asymptotische Näherung der Besselschen Funktion 
J,(e) absolut kleiner als 


st  const |r |"?/2 


2 
ist. Demnach ist der Rest beim Pol == r auch beschränkt. 


f) Das Potential Ö, im Gebiet kleiner Entfernung. Nach der bekannten Be- 
ziehung 


oo 


1 
—h+2)& 7 dE = ———— 0 ....... 0.0. (40) 
ik (Er) 3 (Afajhr 
0 
kann man schreiben 
& 1 w? Jo(Er) 
I= | e%+at 7 ER 2 + III IIE. . . (4). 
a a u a 
g 9 
0 
Setzt man ferner 
2 
san mn? [err-ea 
0 
so folgt 
TerlverWr ine Joe r) „d£—=Re 2 eh+z-ir- sinn 4 dpt. . (42) 
u 7 Pr 
g g 
0 0 


Mit der neuen Abkürzung =h-+2z—i:r-sinp wollen wir die Integration der rechten 
Seite von (42) zuerst nach der Variablen £ ausführen. Den Integrationsweg wollen wir in der 
komplexen £-Ebene so wählen, daß er aus der reellen Achse von Null bis R besteht, wobei 


2 2 
R> = ist, und der singuläre Punkt &= ” mittels eines halben nach oben gerichteten Kreis- 


ER 
* ” * h ” In ” ” 
bogens umgangen wird, ferner aus einem Kreisbogen von Rbis R-e * und einer Geraden, die 
Ei3 


ja 
von R'e * wieder nach Null zurückführt. Da in dem umlaufenen Gebiet der Integrand sich 
überall analytisch verhält, so ist nach dem E Nager Integralsatz das Umlaufintegral 


gleich Null. Auf dem Kreisbogen von Rbis R-e “ bleibt Re(a)> 0, so daß der Betrag == ent- 
sprechenden Integrals mit R— o verschwindet. Wir bekommen also 


el 2? ee a Ser (43) 
ne RER 
. g R g 
;? w? 
Wenn man weiter age re =. 7 =ß und wu—f=t setzt, dann ergibt sich 
cos 27 [+7 oo 
d 
Is —a$ 2. = — BETEN (44). 
Be 
g 
ne % 
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Für den Hauptwert des „„Exponentialintegrals“ auf der rechten Seite von (44) führen wir die 
Bezeichnung ein 


Für die m-te Ableitung der #i-Funktion gilt folgende Reihenentwicklung: 


B.im(f) = m-nte I nr 


2 2 
Wir bilden zweckmäßig die Taylorreihe um die Punkte (h+2)- und En r-sino; denn 
g 9 


für diese Punkte ist Ei(f) aus Tabellen erfaßbar. Somit erhalten wir 


; | ‚wo? \m 
aT Sa\ ee +) „sinn p/ 
m=1n=1 (m—n)! \#+9) = mM 
g 


Ei remis) =Ei a4 — 


Diese Reihe ist für die weitere Integration nach der Veränderlichen 9 ohne Schwierigkeit zu 
gebrauchen. Damit lautet 


a? e oo m N u ei ‘m 
I; ar et g il) | + ar or > en Br n — ? 
re 


und 


a I4m 
De g tsir sn -(i N— Bil(h+2) “+, > > wo r m p (43a). 
g £ 


=0 m=1n=1 


Jetzt wollen wir die Integration in (42) nach der Varibalen y ausführen. Für die Inte- 
2 


gration des halben Residuums des singulären Punktes &= ® in der komplexen &-Ebene und des 


2 
Koeffizienten von Ei a4) “ in (43a) benutzen wir die Reihenentwicklung nach den Bessel- 


schen Funktionen”) 


a 2 = 2 N r 
ne lEr)+2 29 r)-cos@un)+ 25 2, Jans (Sr) -sin(n+) 


und damit erhalten wir 


jo -31& r)+ ler ) Rt (45). 
J | 


»Das letzte Glied von (43a) ist ein Polynom von sin, das für die nahe. Umgebung: des 
Koordinatenursprungs gut konvergiert. Für die Integration von diesem Glied benutzt man 
vorwiegend die beiden folgenden von Wallis eingeführten Beziehungen 


E47 


2 
2:4-6---2n 


a Da u a er, 
IE Be Ag ea 
(n>0, ganz). 


0) 

BR 1:3-5:(@n—1). 
oe. 1.35:"@n—1) # 

[sero Le 2:4:6:-2n 2 


. 12) Webster-Szegö: Partielle Eeeiiunden der mathematischen Physik, 8.373, 1930. 
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Damit können wir schließlich den vollständigen Ausdruck des Potentials &, aufstellen; er gilt 
im ganzen Flüssigkeitsbereich. Mithin wird 


= Q oa 
4 EilzpEz 
ww? er g: 
Ss 2 2 r). Bil+9® + Zmrmuls NE: (46). 
e 5 n]2 
Be ee Ben -- —r > 
+ Re 20° (-=}) ann 9 sin’ tm. dp cos wt 


nn 


g) Die stehenden Wellen im Gebiet kleiner Entfernung. Für die Oberflächen- 
erhehung der Flüssigkeit im Gebiet kleiner Entfernung, die von ©, verursacht wird, erhalten 
wir gemäß Gleichung (3) 


u mtntmsnot a. 


worin 


ne): 


en g T 


REN BET san inltmm » 
m—n)! m(1l) = ) 908 
g 0) 


S 

| 

> 
ale 
Y 
IV. 
I 


Streng genommen gilt Gleichung (47) für den ganzen in Frage kommenden Flüssigkeitsbereich 
Aus Konvergenzgründen aber können wir (47) praktisch nur auf das Gebiet kleiner Entfernung 
von dem Koordinatenursprung anwenden. Infolgedessen wurde Z,, zwecks Unterscheidung von 
&,, mit einem Index % (k bedeutet das Gebiet kleiner Entfernung) versehen. 


Man beachte, daß bei der Auswertung von u, nur gerade Zahlen von [+ min Frage kommen, 
da ungerade Zahlen imaginäre Werte bewirken. Mithin hat man 


v! 


| IR 15 ru ey 1% 1:3:5:7 Say rat 8 


5 2.46% 7 2-.4:.6-.8 
Sinngemäß läßt sich diese Reihe ohne Schwierigkeit für höhere Werte von l hinschreiben. Für - 


2 
den festgelegten Wert 7 h=1 vereinfacht sich die Summation der erhaltenen Reihen: 
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4 6 
6 ) + 0,15625 7) —.0,06682 7) + 0,04001 3) —.0,02773 7 
i 16 
L 0.020087, Y — 0,01612 (5) + 0,01309 | 3) — .0,01089 volEas + 0,00925 05) 


fs = —0,50000 


Se 


24 26 
— 0,00778, a 10 00007, ,) — 0,00618 u 1.0,00543 e ; ne —.0,00488 7)" 
( 


+ 0,00442 GN EHE, 


Die Reihe (48) konvergiert bis zu dem Wert = 1 gut. Die Auswertung von u, und u, ist an 

Hand von Tafeln schnell auszuführen. Wie bei der Untersuchung für das Gebiet großer Ent- 
3 

fernung, setzen wir auch In Nachdem u,, 4, und u; einzeln ausgewertet sind, muß man 


sie noch mit dem gemeinsamen Faktor B= multiplizieren. 


wo 
2ngh 
Fran are Tee 


SSSENEIEMES 
"NEAEDARRER 


Y- 


Bild.6. Momentbilder für die Oberflächenerhebungen £, 7, und Eaoin Abhängigkeit von der Eon T 
für den Zeitpunkt wi = 3 E17 


Das errechnete Ergebnis £,; ist auf der linken Seite des Bildes 6 graphisch aufgetragen. In 
demselben Bild findet man rechts auch einen Teil von £,, (9 bedeutet große Entfernung). Bei- 
spielsweise bei r=10 cm beträgt die Differenz [E. 16 —&, | = 1,396 —1,143= 0,253 cm, und wenn 

sich dieser Unterschied auf Z,, bezieht, beträgt sie etwa 18%. Das heißt, die Näherungsformel 
- für Z,, ist für die Entfernung r="h unter Umständen nicht recht befriedigend. 

Aus (47) ist es auch ersichtlich, daß der Punkt z=0, r=0 nur dann singulär wird, wenn 

die ‚Quelle der Oberfläche unendlich nahe rückt. Im ganzen Gebiet ist u, unabhängig von der 


1 
ni eguenz und Ha, klein gegen Ih. Das Glied ee keinen Wellenchärakter; es el 
En) 


nur die Rolle einer mit der Quelle phasengleichen Grundform der Oberfläche. Die Wellen- 
 erscheinung ist daher hauptsächlich durch die Überlagerung von Besselschen Funktionen 
nullter und verschiedener ungerader Ordnungen charakterisiert. Ferner läßt sich noch fest- 
E | | Bl . (o: ”) $ 3 g 3/2 
Br stellen, daß aus der Se = Pr 7 ie Grenzbeziehung \w°e er. 
ne, lolgt, Wenn der Quellpunkt an die freie Oberfläche rückt, wird Bu, im Gebiet kleiner Ent- 
fernung der dritten Potenz der Kreisfrequenz © proportional sein. 
h) Die fortschreitenden Wellen im Gebiez kleiner Entfernung. Man kann 


E Ei ferner noch zu dem Ausdruck c, (47) die Beziehung 
Be. : e E E aeple os: 


(48). 
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Kleine Mitteilungen 
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hinzufügen, um die fortschreitenden Wellen zu bekommen. Mithin hat man 


el are g? 


ge 2 
RER [300501520 (ntntm) sin mt .".. (49). 


Der Ausdruck (49) läßt sich schwierig zu einer einfacheren Form bringen. Für numerische Aus- 
wertung bietet er aber keine besondere Schwierigkeit. 

Hat man nicht eine einzige kugelige Quelle unter der Oberfläche, sondern stetig längs der 
senkrechten Linie verteilte Quellen, so kann man die Integration in änlicher Weise durch- 
führen. Für ein Quellsystem, das längs einer senkrechten Strecke mit'konstanter Ergiebigkeit 
verteilt ist, hat der Verfasser im Originalexemplar seiner Dissertation (Göttingen 1945) die 
Rechnung durchgeführt. Aber auch noch bei vielen andern Verteilungsgesetzen läßt sich die 
Integration in ähnlicher Weise ausführen, wie sie in dieser Arbeit festgelegt wurde. 


Eingegangen am 29. September 1951. 


KLEINE MITTEILUNGEN 


Brückenlösungen beim Problem der achsen- 
symmetrischen Torsion. 


Das Verfahren dererzeugenden Funktionen wurde im 
Falle der torsionsfreien achsensymmetrischen Defor- 
mation wiederholt erfolgreich angewendet; u. a. wurde 
vom Verf. für das in Frage kommende Bipotentialpro- 
blem bei Verwendung von Polarkoordinaten ein Sy- 
stem von Brückenlösungen entwickelt 1), wie sie im 
Fallederinnen und außen belasteten Hohlkugel heran- 
zuziehen sind. Die vorliegende Arbeit richtet sich 
auf die analoge Aufgabe beim Torsionsproblem und 
gibt damit Ergänzungen zu älteren Arbeiten von 
A.Föppl?), Fr. A. Willers?) und des Verf. ?). 

Bei fehlender Massenkraft haben wir unter Anwen- 
dung von Polarkoordinaten o, x, für den Fall, daß 
alle Spannungen und Verschiebungen von der Meri- 
diankoordinate ® unabhängig sind, als einzige den Tor- 
sionsfall charakterisierende Gleichgewichtsbedingung?): 


9% 1 07x 1 
- +- 2 +-—-|3%-+27ry cot =0... (1). 
90 th e +27 8x) (1) 
Darin bedeuten ro, 7 die Komponenten der Schub- 
spannung, die im Meridianschnitt in der Richtung der 
wachsenden o- bzw. x-Werte wirken. Für das Ober- 
flächenelement des Umdrehungskörpers mit der nach 
außen gerichteten Normalen » ergibt sich als Dreh- 
spannung: 
tH= 70608 (0,9) + 7208 (4,9). » 2... (2) 
Bezeichnet v die lineare Drehverschiebung irgend- 
eines Punktes des Rotatlionskörpers, so lauten die Be- 
ziehungen zwischen Spannung und Formänderung ®): 


= G@osin RERLREAN = (sin AN: 
ig 77 eainy) TAX, esinyx 
(3) 
Die Grundgleichung des Torsionsproblems ergibt sich 
durch Einsetzen der Ausdrücke (3) in die Gleich e- 
wichtsbedingung (1). Sie nimmt bei Einführung der 
Drehungskomponenten 


1 0) 
200 = —— —— (bosiny); 
eT pesiny 0, | ee) (4) 
u nom | 
*  osiny de | Page 
die Gestalt an; | 
PINS ne) 


0% do 


i RT ER. 30 (1950), 8. 50 ff. 
nchen Akad. d. Wiss. Sitzungsber. 35 (1905), 8. 249, 504; 
Zeitschr.d. V.d. I. 50 (1906).8. 1082. N GR 
#) Zeitschr. Math. Phys. 55 (1907), 8. 225—263. 
*) Z. angew, Math. Mech. 11 (1931), 8. 8—15. 
°)A.E.H. Love, Elastizität 1907, 8. 108. 
°)A.E.H. Love, Elastizität 1907, 8. 67. 


wird also durch den Ansatz befriedigt 


BE Ar 2 ‚ 
Nr ö 2x Late A 
Koppelt man diesen Ansatz mit (4), so ergibt sich: 
0 Ser ee 
7, Pesind=e sın X do ’ 
ö 2 7% OR 
Zr a N 3 


woraus i 
DT N OT 
DE (e sin x tin) = 5 


die erzeugende Funktion p genügt also der 
Laplaceschen Gleichung’). 

Nun gibt es zu jeder einem achsensymmetrischen 
Problem entsprechenden Potentialfunktion p eine 
assoziierte „Strömungsfunktion" y, die mit ihr durch 
die Beziehungen 


ler en: 5 
verknüpft ist. Gemäß (4), (4) ist im vorliegenden Fall 


y = vosiny. Ferner kann man?) ein derartiges Funk- 
tionpaar sich verschaffen durch den Ansatz: 


0 i 
o--"5, v=esing- da, ..:.(b) 


wo 9, eine achsensymmetrische Potentialfunktion. 
Im Falle des Problems der achsensymmetrischen Tor- 
sion gewinnt man dann anschließend®) die Drehver- 


schiebung vd = 2 — ferner v = BASE und die 
osin x 


osiny; 
Spannungskomponenten 
u ov 
%o=@ INN, = j =, 
f) esinx d4 %=@siny 2% 


Wie beider torsionsfreien achsensymmetrischen 
Deformation von Umdrehungskörpern beanspruchen 
besondere Aufmerksamkeit auch die für den Torsions- 
fallzu gewärtigenden Brückenlösungen, bei denen 
Potenzen von g mit den um nullliegenden Exponenten 
als Faktoren auftreten. Natürlich scheiden hier die 
nur der Bipotentialgleichung, aber nicht der Poten- 
tialgleichung genügenden Funktionen aus. Das sind, 
wenn wir die Lösungsbezifferung benutzen, die vom . 
Verf. in dem unter!) zitierten Aufsatz angewendet 
wurde, die folgenden: 


4,, B,, B;, D, D, H,, H,, H, T,, U,, U. 


?) Vgl.}). 
®) Vgl.®). 


© u, 
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Es verbleiben somit die nachstehenden Lösungen: 1 1 ; 
[ro = a sin, 


cos 0 
[Bılpı= 5; le ö 1 
eg e 4 e 9= no muy=-—--; 
ende ng e 8 
al X du Dr rn, = — co8X 
i 1 1 cosyx 
; y—-——cegy; we——— 3 
er een De EN Kr BR 
osinx e esiny @ Wa Urkoeg 
. 0 =. 6 . si — 
%=Gesiny >= —30 8; 0? } 0? sin?y 
5 [U] 9 =Ino + Insiny; 
= Gsinyn—0. = —Ino—1-—Insiny v=ocose 
* b=ctgy; ER T Fe ; 
Diese Lösung wurde zuerst von A. Föppl°) auf er ER oz 


anderm Wege abgeleitet. 1 —sin®y —2sinx cos?y 


1 u=Gsiny: - 
[D,] m inted; e-- [mis)=0; i “ sint y 
: RT Eo8% 
ea hie Er 2 ep > sin? \ 
y=osiny Den y 15 Ir X ; 
Inversionvon Torsionszuständen. 
= - ; = ee ; Wie beim ebenen Spannungsproblem gibt es auch *) 
esıny 0? sin? y bei der achsensymmetrischen Torsion die Möglichkeit, 
2G 260 cosx eine neue Spannungsverteilung durch Inversion abzu- 
ER ee a er ze leiten. Bekanntlich erhält man durch die Substitution 
0? sinx 0? sin?x 


91 = 09, aus einer harmonischen Funktion 9, eine 
neue Funktion 9], die wiederum der Laplaceschen 
[23] 9 =Intg 93 Gleichung genügt, wenn man gleichzeitig die Punkte 


des Raums durch reziproke Radien abbildet: 0’ =, 5 


X=yx; 6 =. Einer Strömungsfunktion 


0) 
1 1 a 001 
b=siny’ ”” osiny’ DEE 
G 20 cosx ordnet sich bei diesem Prozeß die Funktion 
mn HOT my‘ a N BE 
x 2 “ 0% % 0 
HA. = ft -I I zu, die im invertierten Raum, also rücksichtlich der 
en 0° (oo, Rense 2 & ; Koordinaten o, x, % wieder die Bedeutung einer Strö- 


mungsfunktion hat. Im invertierten Raum gilt daher 
folgendes, einen Torsionszustand charakterisieren- 


1 % 
= — Xlntg- +1); 
“ o* (or Ku ) des Formelsystem: 
_sinx(_ . r2 : d 
re | sin gm tg &+etez); v--0,(8); = 0%; 
ya — sinyIntg +ctgy); EEE AEN al». 
® 2 ST 2 sin? y/' 
rk Ba ers \, 1. Bsp. Aus Lösung H, ergibt sich durch Inversion 
TE 0’ (im as 2 = sin? :) £ folgender Torsionstfall: : 
3@/. , 
a  (inzimts$ +0tgz):; - 1-30 jmemteh tanz]; 
2G | y r SE) 49.6084 ]° 
Ei er % =20@ctg!y + 34, |- Inte. + | 
IT) a=meo+in(l+cosy); 2. Bsp. Aus Lösung T' folgt durch Inversion der Tor- 
e=—-Inoe—1-—In(1+cosx); sionsspannungszustand: 
_ wi z aina d . ! 2% v4 
Yy=TeTLcosg’ n=Goten; 
i 1 
2 RE UL 2 oe 
a 1+cosx ’ 3 e(l-+cos x) ’ = + 30m (das; 
| Be Se er \ q IR fe} * —eosiny. 
WO Tr con) 0 9° = Go’tg a3 GSX TH cosy)® 
fe 
u=- 18 5- x > Get -300t0 5. 
” Sitzungsberichte der math., phys. Klasse der bayr. Ak. 35 (1905) Berlin-Schlachtensee. A. Timpe. 


249 u.S8. 
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BUCHBESPRECHUNGEN 


Herbert Goldstein (Ph.D., Harvard University), 
Classical Mechanics. 3098. Cambridge 
(Mass.) 1951. Addison-Wesley Press, Inc., Preis 7,—$. 

Mit dem vorliegenden Buche wendet sich der Verf. 
vorwiegend an Studierende der theoretischen Physik 
und Mathematik, wobei der neueste Stand dieses 
Fachgebietes noch nicht Berücksichtigung findet, 
sondern der Stoff im wesentlichen nur die klassischen 
Bestandteile der theoretischen Mechanik umfaßt. 
Andererseits gewinnt der Leser ein außerordentlich 
umfangreiches und gut durchgearbeitetes Rüstzeug, 
das bei der Beschäftigung mit schwierigeren Teilge- 
bieten der modernen theoretischen Physik sehr vor- 
teilhaft ist. Das Buch gliedert sich in folgende Ab- 
schnitte: Elementare Grundlagen, Variationsprinzi- 
pien und Lagrange sche Gleichungen, Zwei-Körper- 
Problem, Kinematik des starren Körpers, Dynamik 
des starren Körpers, spezielle relativistische Mechanik, 
dieHamilton schen Gleichungen, kanonische Trans- 
formationen, die Hamilton-Jacobische Theorie, 
kleine Schwingungen, Lagrange sche undHamilton- 
sche Formulierungen für Kontinua und Felder. Die 
Darstellung läßt erkennen, daß der Verf. den Stoffvor- 
wiegend von der mathematisch-formalistischen Seite 
her behandelt, während die physikalische und tech- 
nische Anwendung mehr im Hintergrund bleibt. So 
sind auch nur wenige Abbildungen beigefügt. Dennoch 
sind eine Reihe von Übungsaufgaben angegeben, teil- 
weise auch mit ausführlichen Erläuterungen, sodaß 
auch der noch ungeübte Leser die Möglichkeit erhält, 
einzelne Probleme selbständig zu bearbeiten. An 
vielen Einzelheiten erkennt man, daß der Verf. be- 
müht gewesen ist, das jeweilige Problem so ausführlich 
zu formulieren, daß die Anwendung für besondere 


Fälle der modernen Physik keine Schwierigkeiten mehr 
bedeutet. Dies ist z. B. bei der Zentralkraftbewegung 
bei der Hamilton-Jacobischen Theorie und bei 
der speziellen Relativitätstheorie der Fall. Das Buch 
kann daher den Studenten und Physikern als zuver- 
lässiges Übungs- und Nachschlagewerk warm emp- 
fohlen werden. 
Dresden. { H.Neuber. 
Dr.-Ing. habil. H. Müller, Mechanik für In- 
genieure. Bd.3:Festigkeits- und Ela- 
stizitätslehre, 254 S. m. 268Abb. Braun- 
schweig-Berlin-Hamburg-Kiel 1951. Georg Wester- 
mann-Verlag. Preis geh. 16,60 DM., geb. 18,60 DM. 


Das vorliegende Buch gehört zur Gruppe der nach 
dem Kriege in größerer Zahl erschienenen Elementar- 
Lehrbücher, welche der regen Nachfrage entgegen 
kamen, die durch die Kriegsverhältnisse und den ein- 
getretenen Mangel an Schrifttum entstanden war. 
Inzwischen ist aber bereits eine Reihe wertvoller Lehr- 
bücher auf dem Gebiete der Mechanik erschienen. Im- 
merhin zeichnet sich das Buch durch sorgfältige Dar- 
stellung und geschickt ausgeführte Abbildungen aus 
und kann für den Anfänger als gute Stütze dienen. 
Der Inhalt umfaßt folgende Teilgebiete: Grundbe- 
griffe, einfache Beanspruchungen, zusammengesetzte 
Beanspruchung, Formänderungen und Formände- 
rungsarbeiten, statisch unbestimmte Systeme. Das 
Buch enthält zahlretche einfache Beispiele, die den 
Leser zum selbständigen Nachrechnen anleiten, wobei 
sich der Verf. jedoch nur wenig an praktisch vor- 
kommende Fälle anlehnt. 


Dresden. H.Neuber. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


NACHRICHTEN 


Zu Prof. Schulers 70. Geburtstag 


Am 5.Febr. dieses Jahres beging Prof. Dr.-Ing. 
Max Schuler, Göttingen, in bester körper- 
licher und geistiger Frische seinen 70. Geburtstag. 
Als Prof. Schuler 1947 emeritiert wurde, konnte 
er. auf eine langjährige erfolgreiche Tätigkeit als 
akademischer Lehrer an der Universität Göttingen 
und als Leiter des seit 1940 ausgebauten neuen In- 
stituts für angewandte Mechanik zurückblicken, das 
nachher bekanntlich den dortigen physikalischen 
Instituten angegliedert‘ worden ist. Der Fachwelt 
und insbesondere seinen Schülern ist Prof. Schuler 
ein Begriff durch die seltene Verknüpfung von 
Wissenschaft und technisch-mechanischer Erfindungs- 
kraft, die seiner Arbeit das Gepräge aufdrückt und 
ihn zweimal zu außergewöhnlichen Erfolgen geführt 
hat. - 

In der ersten Phase seiner Laufbahn ging es um 
die Anwendung des Kreisels zu Zwecken der Navi- 
gation und Steuerung. Nach abgeschlossenem Stu- 
dium an der T. H. München betätigte er sich seit 1907 


“als Mitarbeiter von Anschütz -Kiel bei der Kon- 


struktion des Kreiselkompasses. Er führte hierbei das 
Prinzip des 84-Minutenpendels ein, konstruierte 1910 
den Einkreiselkompaß, klärte den Schlingerfehler bei 
Kreiselkompassen auf und schuf schließlich 1912 den 
verbesserten ‚Dreikreiselkompaß‘“. Weitere Erfin- 
dungen Schulers aus dieser Zeit sind der Flieger- 
horizont und der Wendezeiger für den Blindflug. 


Herausgeber und Hauptschriftleiter: Prof. Dr. Fr. A. Willers, Dresden. Verlag: Akademie-Verlag GmbH., Berlin NW 7, Schiffbauer- 


Die zweite Phase der Laufbahn Schulers wurde 
durch seine Habilitierung für angewandte Mechanik 
in Göttingen 1924 eingeleitet. In der Folgezeit wandte 
er sich dort neben seiner Lehrtätigkeit in Zusammen- 
arbeit mit der Sternwarte dem Problem der Präsizions- 
messung von Zeit und Schwere zu. Beginnend mit 
dem Ausgleichspendel (1929) gelang ihm nochmals in 
stiller und zäher Arbeit eine Entwicklung, die wie das 
frühere Beispiel des Kreiselkompasses zeigt, daß die 
Mechanik auch heute noch zur Zusammenarbeit und 
zum Wettstreit mit den moderneren Zweigen der 
Physik fähig ist. Heute besitzt die Wissenschaft in der 
Schuleruhr ein Hilfsmittel, welches nicht nur die 
genauesten Zeitmessungen ermöglicht, sondern im 
Verein mit der Quarzuhr es auch gestattet, die 
anderen beiden Größen der Pendelgleichung, nämlich 
Fallbeschleunigung und Pendellänge der genauesten 
Kontrolle zu unterziehen. Dadurch ist einerseits der 
Beitrag der Gravitation des Mondes zur Pendelbe- 
schleunigung zutage getreten und sogar die elastische 
Verformung der Erde durch die Anziehung des Mondes 


in den ge a gerückt, andererseits - 


beobachtete Schuler sprunghafte Längenände- 
rungen des Invarstahls für die Pendelstange, die 
wichtige Folgerungen für das molekulare Verhalten 
der Invarstähle und anderer Werkstoffe und für deren 
Verwendung im Vermessungs- und Markscheidewesen 
erwarten lassen. 


Landshut. H. Gebelein. | 
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